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•  /  when  applied  to  explosions  in  other  regions  of  the  United  States  and  in 
the  French  Sahara,  gives  realistic  yield  estimates.  The  technique  also 
is  applied  t^- selected  Soviet  explosions  in  East  Kazakh.  m^CLg)  and 
tBjj(P)  values  were  used  to  estimate  the  mtj(P)  bias  between  NTS  and 
eastem^North  America.  Assuming  that  explosions  and  earthquakes  of  the 
same  mb(P)  value  excite  Lg  waves  of  equal  amplitude,  the  P-wave 
magnitude  bias  between  NTS  and  eastern  North  America  is  0.31  magnitude 
units.  A  tentative  value  for  the  bias  between  NTS  and  Shagan  River  is 
0.41  magnitude  units,  but  this  value  may  be  changed.. when  further 
study  of  anelastic  attenuation  of  Lg  in  Eurasia  is  completed. 

Tbe  effects  of  velocity  structure,  crustal  Q,  and  sediment  Q  values 
on  the  attenuation  of  Lg  waves  were  studied  using  theoretical  amplitude 
spectra  and  synthetic  seismograms  for  various  models.  Comparisons  with 
observations  suggest  that  crustal  velocity  structure  can  have  a 
significant  effect  on  Lg  attenuation  and  that  low-Q  sediments  have 
a  profound  effect  of  the  variation  of  Lg  attenuation  in  stable  regions. 
By  contrast,  in  areas  where  crustal  Q  values  are  low,  there  is  little 
effect  of  sediment  Q  values  on  the  attenuation  of  Lg.  Low  Q  values 
for  the  crust  of  the  western  United  States  seem  to  be  required  in  order 
to  match  the  variation  of  frequency  in  the  wave  forms  of  Lg. 


Frequency-dependence  of  crustal  Q  seems  significant  in  regions  of 
high  Q,  but  are  small  or  non-existent  in  regions  of  low  Q  values. 

/f etfu'oxJs  A;-i  '  dt,  ^  St>r  \  m, 

etpycS’/c  QwJ  Seismic  //€/d. 


5 


]e<ju'tn< 
c 


J) 


W it  z-Lp  r 


(I 


Accession  For 

GRA&I 


n:  :  T4B 

jjiced 

t.,,u f  :<  fit.  Von. 


i-"  ; at  ion/ 

;■  v ‘.  ' -t-  Itv  Codes  j 
v  .  :  or..l/or 

Dlst  •  v  ■- Val 


(S5. ) 


— \ 


SECURITY  CLASSIFICATION  OF  THIS  PACEmwi  DM*  Enfnd) 


i 


AFOSR-TR.  8  4-1 

Final  Technical  Report 
1  October  1982  -  30  September  1984 


ARP  A  Order: 

Program  Code: 

Contractor : 

Effective  Date  of  Contract: 
Contract  Expiration  Date: 
Amount  of  Contract: 

Contract  Number: 

Principal  Investigators: 

Telephone : 

Short  Title: 


4397 
3D  60 

Saint  Louis  University 

1  October  1982 

30  September  1984 

$279,034.00 

F49620-83-C-0015 

Brian  J.  Mitchell 
Otto  W.  Nuttli 

314-658-3131 

Surface  Wave  Attenuation 


Sponsored  by 

Advanced  Research  Projects  Agency  (DOD) 

ARPA  Order  No.  4397 

Monitored  by  AFOSR  under  Contract  //F49620-83-C-0015 


The  views  and  conclusions  contained  in  this  document  are  those  of 
the  authors  and  should  not  be  interpreted  as  necessarily  representing 
the  official  policies,  either  expressed  or  implied,  of  the  Defense 
Advanced  Research  Projects  Agency  of  the  U.S.  Government. 


Approved  for  public  release 
distribution  unlimited. 


Table  of  Concents 


Technical  Report  Summary 

Methodology  of  using  Lg  Waves  for  Estimating 
Explosion  Yield  and  Body-Wave  Magnitude  Bias 
between  Test  Sites 

Coda-Q  Studies  for  the  Indian  Subcontinent 

Attenuation  of  High-Frequency  Earthquake  Waves 
in  South  America 

Effect  of  Crustal  Velocity  Structure  and  Crustal 
Q  Values  on  the  Amplitudes  and  Wave  Forma  of  Lg 

Summary  of  Studies  on  the  Frequency  Dependence 
of  Qg  in  the  Continental  Crust 


*0‘ 

ItA* 


TECHNICAL  SUMMARY 


Lg  Yield  And.  Magnitude  Bias  Esfr-liBafcaa 

The  methodology  for  using  1-see  period  Lg  amplitudes  has  been 
refined  for  determining  yields  of  underground  nuclear  explosions  at  NTS. 
In  particular,  special  attention  vas  devoted  to  obtaining  the  best  pos¬ 
sible  estimates  of  the  coefficient  of  anelastic  attenuation  for  paths 
from  NTS  to  the  WWSSN  stations  BKS.  DUG,  and  TUC.  In  the  first  itera¬ 
tion  the  coda-Q  method  is  applied  to  determine  Q  (1-Hz  value  of  Q)  and 

o 

the  parameter  $  which  is  a  measure  of  the  frequency  dependence  of  Q. 

Then  small  corrections  to  the  Q  values  are  made  so  as  to  minimize  the 

o 

deviations  between  the  station-average  m^Lg)  and  the  individual  m^Lg) 
values.  Two  problems  were  encountered  and  overcome  in  applying  the 
coda-Q  method.  The  first  concerns  the  interference  of  low-frequency 
fundamental-mode  Rayleigh  waves  with  the  Lg  coda  waves  about  1  to  3 
minutes  after  the  onset  of  Lg,  which  if  not  corrected  for  will  give  fre¬ 
quency  values  that  are  too  low  at  relatively  small  lapse  times.  The 
second  concerns  the  leveling  off  of  frequency  for  large  lapse  times, 
with  the  lapse  time  at  which  the  frequency  attains  a  constant  value 
being  a  function  of  the  magnitude  of  the  explosion.  This  phenomenon  can 
Incorrectly  result  in  too  large  estimates  of  wave  frequency  at  long 
lapse  times.  Both  phenonmena,  especially  the  latter,  can  result  in 
underestimates  of  $  and  overestimates  of  Qq  values.  Nevada  earthquakes, 
for  which  both  m^(P)  and  m^(Lg)  were  obtained,  support  the  Qq  and  Jj 
values  found  in  this  study. 

Lg  amplitudes  and  m^(Lg)  values  were  obtained  for  all  NTS  explo¬ 
sions  of  announced  yield,  and  calibration  curves  are  presented  for  shots 


in  hard  rook  aad  in  alluvium.  The  standard  deviation  of  the  hard  rock 
m^(Lg)  values  is  0.06  magnitude  units  (equivalent  to  about  a  15$  uncer¬ 
tainty  in  yield  value). 

Lg  amplitudes  also  were  used  to  estimate  yields  of  non- NTS  events 
in  the  United  States  (SALMON,  GASB00GT,  RULISON  and  RIO  BLANCA)  and  in 
the  French  Sahara  (SAPHIR  and  RUB IS) .  The  hard-rock  NTS  calibration 
curve,  vhen  applied  to  these  data,  gave  yield  values  close  to  the 
announced  one.  This  suggests  that  there  is  no  magnitude  bias  in  m^(Lg), 
in  contrast  to  observed  bias  for  m^P),  and  that  the  excitation  of  Lg 
waves  for  an  explosion  of  a  given  yield  is  the  same  in  all  continental 
areas. 


Q0«  (j  and  m^(Lg)  values  were  obtained  for  Soviet  explosions  at  the 
Shagan  River  and  Degelen  Mountain  regions  of  the  East  Kazakh  test  site. 
Preliminary  results  indicate  a  yield  of  200  kt  for  the  largest  explosion 
considered,  that  of  Septemver  14,  1980.  For  the  only  Soviet  explosion 
of  announced  yield  at  East  Kazakh,  that  of  January  15,  1965,  the 
announced  yield  was  125  kt  and  the  m^(Lg) -estimated  yield  was  105  kt. 

The  mb(P)  and 

m^P)  bias  of  NTS  relative  to  eastern  North  America  is  0.31  magnitude 
units  or  larger,  depending  on  whether  the  excitation  of  Lg  waves  is  the 
same  for  explosions  as  for  earthquakes  of  the  same  m^(P),  or  whether  the 
Lg  excitation  for  explosions  is  less. 

Tentatively  the  m^(P)  bias  between  NTS  and  Shagan  River  events  is 
found  to  be  0.41  magnitude  units.  However,  the  Qq  and  $  values  for  the 


m^(Lg)  values  for  NTS  explosions  indicate  that  the 


Soviet  explosions  need  to  be  further  studied,  which  will  be  done  in  the 


work  to  follow.  If  the  Shagan  River  P-wave  attenuation  in  the  astheno- 
sphere  is  the  same  as  for  eastern  North  America,  the  bias  between  it  and 
NTS  will  be  0.31  magnitude  units. 


Lg  Attenuation  In  South  America 

The  coda-Q  method  was  applied  to  earthquakes  occurring  in  South 
America  and  recorded  by  WWSSN  stations  across  the  continent.  Amplitude 
attenuation  was  used  as  an  independent  method  of  obtaining  the  same 
information. 

Both  methods  give  low  Q  values  (135  to  320)  in  the  northern  and 

o 

wstern  coastal  areas  of  South  America,  the  regions  of  neotectonics.  The 

shield  and  geologically  older  areas  of  eastern  South  America  have  an 

average  (3  value  of  about  700.  Regions  of  low  Q  have  6  values  of  0.4 
o  O’ 

to  0.7 •  and  regions  of  high  Qq  have  $  values  of  0.0  to  0.2. 

Lg  Attenuation  la  India  and.  Pakistan 

The  coda-Q  method  was  applied  to  earthquakes  occurring  in  India  and 
recorded  by  WVSSN  stations  in  India  and  Pakistan.  Amplitude  attenuation 
was  used  to  confirm  the  coda-Q  results  for  selected  stations. 

The  mountainous  areas  of  northern  India  and  Pakistan  show  rela¬ 
tively  low  Qq  (300)  and  relatively  high  ^  (0.4)  average  values.  The 

shield  area  of  northern  India  shows  relatively  high  Q  (700)  and  rela- 

o 

tlvely  low  $  (0.2)  values. 

The  research  of  B.J.  Mitchell  has  been  concentrated  on  two  major 
aspects  of  seismic  wave  attenuation.  First,  several  factors  which 
influence  the  attenuation  of  seismic  waves  at  regional  distances  have 


I%»  »I.%»I  |  %« I  ,7m  InlV 


been  thoroughly  Investigated.  These  Include  crustal  velocity  structure, 
crustal  Q  structure,  and  Q  values  and  thicknesses  of  surfioial  sedi¬ 
ments.  Second,  the  frequency  dependence  of  shear  wave  Q  in  the  crust  of 
South  America  and  India  have  been  studied  using  both  1-Hz  Lg  and  Lg 
waves  at  Intermediate  periods.  The  results  of  these  studies  can  be  sum¬ 
marized  as  follows. 

£f£e££&  at  Crustal  Velocity  Structure  fin  Lg  Attenuation 

In  the  eastern  United  States  the  presence  of  low-velocity  sediments 
contributes  to  the  length  of  the  Lg  coda  and  smooths  the  wave  form  which 
for  simple  crustal  models  of  few  layers  and  no  sediments  is  splkey  in 
appearance.  Known  differences  in  velocity  structure  between  the  eastern 
United  States  and  Basin- and-Range  province  can  lead  to  differences  in 
mb(Lg)  predicted  from  synthetic  seismograms  computed  for  those  models. 
Assuming  identical  Q  models  for  the  two  regions,  we  find  magnitude 
differences  of  0.2  to  0.3  units  for  1-Hz  Lg  waves.  Magnitude  differ¬ 
ences  remain  substantial  to  periods  as  large  as  6s,  but  become  much 
smaller  at  longer  periods. 

Effect  fi t  lfiVrQ  &edJjnej&a.  fin  Lg  Attenuation 

The  large  variation  observed  for  Lg  Q  values  at  1-Hz  across  the 
eastern  United  States  can  be  explained  as  being  due  to  variations  in 
thickness  of  low-Q  sediments.  The  frequency-dependent  Q  model  of 
Mitchell  (I960)  covered  with  varying  thicknesses  of  sediments  having 
realistic  values  of  Q  produced  reasonable  agreement  with  observed  val- 
lues  of  Lg  Q  in  Nebraska.  Pennsylvania,  and  gulf  coast  Louisiana,  where 
they  range  between  about  1300  and  400.  In  regions,  such  as  the  western 


United  States,  where  crustal  Q  values  are  thought  to  be  much  lower  than 
those  in  stable  regions,  the  presence  of  low-Q  sediments  has  a  much 
smaller  effect  on  the  attenuation  of  1-Hz  Lg. 

Effect  of  Crustal  £  Structure  on  Lg  Attenuation 

Both  models  with  a  thick  layer  of  low-Q  sediments  overlying  a 
high-Q  upper  crust  and  models  with  a  low-Q  upper  crust  can  explain  the 
attenuation  of  Lg  waves  at  intermediate  periods  in  the  western  United 
States.  The  observed  attenuation  of  1-Hz  Lg  waves,  however,  appears  to 
require  a  model  for  the  Basin-and-Range  which  has  low  Q  values  in  the 
upper  crust. 

Frequency-Dependence  of  Apparent  Shear-Wave  £  In  ihe  crust 

In  both  South  America  and  India,  the  upper  crust  in  the  stable  por¬ 
tions  of  those  regions  is  characterized  by  frequency-dependent  values  of 
shear- wave  Q.  In  order  to  model  1-Hz  Lg  waves  for  those  stable  regions, 
shear  wave  Q  values  in  the  upper  crust  must  be  higher  than  those  needed 
to  explain  shear-wave  Q  values  at  intermediate  periods.  By  contrast, 
the  upper  crust  in  the  tectonically  active  protions  of  both  continents 
appear  to  have  low  values  of  shear-wave  Q  which  do  not  depend  upon  fre¬ 


quency 


METHODOLOGY  OP  USING  Lg  WAVES  FOR  ESTIMATING  EXPLOSION 
YIELD  AND  BODY-WAVE  MAGNITUDE  BIAS  BETWEEN  TEST  SITES 

Otto  W.  Nuttli 

INTRODUCTION 

Lg  waves  are  short-period,  higher-mode  surface  waves  that  propagate 
in  the  continental  crust.  They  were  first  described  by  Press  and  Ewing 
(1952),  who  noted  that  even  a  small  segment  of  oceanic  crustal  path  (as 
little  as  200  km)  was  sufficient  to  extinguish  them.  Their  group  velocity 
at  a  period  of  1  sec  is  approximately  3*5  km/sec.  They  exhibit  vertical 
and  radial  as  well  as  transverse  particle  motion,  suggesting  that  they  are 
a  combination  of  higher-mode  Raleigh  and  Love  waves.  Nuttli  (19?3)» 
making  use  of  the  fact  that  Lg  waves  in  general  have  larger  amplitudes  than 
P  waves  at  regional  distances  in  eastern  North  America,  employed  them  to 
estimate  body-wave  magnitude  (m,  )  of  small  explosions  and  earthquakes  for 
which  no  P-wave  data  could  be  obtained.  Furthermore,  their  amplitudes  appear 
to  be  less  sensitive  to  focal  mechanism  and  focal  depth  than  those  of  P  waves, 
which  results  in  a  smaller  uncertainty  and  standard  deviation  in  the  estimated 
m^  value  obtained  from  Lg  amplitudes. 

Lg  ATTENUATION 

In  the  time  domain  at  regional  distances  the  attenuation  of  surface-wave 
amplitudes  due  to  geometric  spreading  and  wave  dispersion  was  shown  by  Ewing 
et  al  (1957)  to  vary  as  &  1  for  normally  dispersed  waves  and  A  for  Airy 
waves,  where  &  is  epicentral  distance.  Observational  data  for  Lg  (Nuttli, 
1973)  suggested  that  the  &  reiation  Is  more  appropriate.  Recently  the 
A  relation  was  confirmed  by  numerical  modelling  of  regional  waves  by 
Campillo  et  al  (1984),  who  showed  that  the  Lg  amplitudes  on  synthetic  seismo¬ 
grams  decayed  as  A  when  no  anelastic  attenuation  (infinite  Q)  was 


assumed. 


Prom  earthquake  data  Nuttli  (1973)  found  that  the  anelastic  attenuation 
of  1-sec  period  Lg  waves  in  eastern  North  America  is  low,  approximately 
0.0006  km~*  (Q  =  1500)9  compared  to  a  value  of  about  0.005  km  *  (Q  =  180) 
for  southern  California.  This  accounts  for  the  fact  that  Lg  waves  are  much 
more  prominent  on  regional-distance  seismograms  for  eastern  North  American 
earthquakes  than  for  western  ones. 

Singh  and  Herrmann  (1983)  used  Herrmann's  (1980)  extension  of  Aki  and 
Chouet's  (1975)  method  for  estimating  Q  from  the  coda  of  Lg  waves  to  map  the 
apparent  Q  of  1-sec  period  Lg  waves  over  the  United  States.  Their  map  shows 
extreme  values  of  140  (California)  and  1350  (central  United  States),  along 
with  considerable  character  or  geographical  variability.  It  suggests  that 
if  Lg  amplitudes  are  to  be  used  to  estimate  explosion  yield  to  an  accuracy 
of  3056  (m^  to  an  accuracy  of  about  0.1),  it  is  inadequate  to  use  an  average 
value  cf  apparent  Q  of  Lg  over  a  region  as  large  as  the  western  United  States. 
Rather  it  is  necessary  to  estimate,  as  accurately  as  possible,  the  apparent 
Q  value  for  particular  source- to- station  paths.  This  statement  takes  on  in¬ 
creasing  significance  for  regions  of  low  Q  (e.g.  NTS)  and/or  for  large  epi- 
central  distances  (e.g.  Soviet  explosions  recorded  by  WWSSN  stations). 
Therefore,  before  Lg  amplitudes  can  be  used  for  determining  m^(Lg)  and  for 
estimating  explosion  ■>*',eld  of  NTS  and  Soviet  explosions,  it  is  necessary  to 
obtain  Q  values  for  specific  source- to- station  paths.  The  problem’s  diffi¬ 
culty  is  compound edby  the  facts  that  the  Lg  waves  do  not  always  have  a 
period  of  1.0  sec  and  the  Q  value  is  frequency  or  period  dependent.  In  this 
report  it  is  assumed  that  the  frequency  dependence  can  be  written  as 
(Mitchell,  1980) 

Q(f )  =  Q0  f  T 

where  Qq  is  the  1-Hz  value,  f  is  wave  frequency  and  J  is  a  number  between 


0  and  1,  given  to  one  significant  figure. 

For  the  purposes  of  this  study  a  first  approximation  of  the  Qq  and  $ 
values  was  obtained  by  using  the  method  of  Herrmann  (1980)  for  individual 
source- to- station  paths,  and  then  by  adjusting  the  Qq  values  ( J  is  not 
allowed  to  vary)  so  as  to  minimize  the  deviations  of  individual  station  m^ 
values  from  the  station-average  value.  Typically  these  adjustments  in 
values  are  found  to  be  less  than  10#  of  the  values  of  the  first  approximations. 

Because  the  m^(Lg)  and  explosion  yield  estimates  depend  critically  on  the 
Q  values  employed,  careful  consideration  must  be  given  to  the  determination  of 
Q.  For  NTS  events,  with  low  Q  paths,  the  WWSSN  station  TUG  as  a  distance  of 
about  700  Jem  was  found  from  observational  data  to  be  the  most  distant  WWSSN 
station  that  gave  consistent  m^(lg)  values.  Data  from  BKS  and  DUG  also  were 
found  to  be  usable.  It  would  be  desirable  to  have  data  from  more  than  three 
stations  for  estimating  an  average  rn^Lg),  but  the  uncertainties  in  amplitude 
associated  with  the  attenuation  corrections  for  longer  paths  outweigh  the 
advantages  of  using  more  of  the  existing  WWSSN  stations.  Therefore,  because 
of  the  limited  number  of  stations,  the  results  of  this  study  are  primarily 
useful  for  developing  a  methodology  of  estimating  m^(Lg)  and  explosion  yield, 
rather  than  of  obtaining  precise  values  of  these  quantities  for  NTS  events. 

In  general  the  Q  values  for  source- to- station  paths  of  Lg  waves  from  ex¬ 
plosions  at  the  East  Kazakh  test  sites  recorded  by  WWSSN  stations  in  Scandi¬ 
navia  and  southern  Asia  are  larger  than  for  the  NTS  region.  Therefore  stations 
at  larger  epicentral  distances  can  be  used  (more  than  4000  km  for  some  of  the 
Scandinavian  stations  for  large  explosions).  But  the  distribution  of  WWSSN 
stations  with  respect  to  the  East  Kazakh  sites  is  far  from  ideal  for  accurate 
evaluation  of  rn^Lg)  and  yield.  Much  better  results  could  be  obtained  from 


seismographs  located  in  the  geological  shield  and  old  plateau  regions  of 
the  Soviet  Union,  at  distances  of  no  more  than  1000  to  2000  km. 


The  Aki  and  Chouet  (1 975)  and  Herrmann  (1980)  studies,  used  in  the 
present  analysis,  consider  the  Ig  coda  to  consist  of  scattered  waves.  There¬ 
fore,  they  can  be  expected  to  give  a  Q  value  that  is  representative  of  the 
region  over  which  they  are  scattered,  rather  than  an  average  value  over  a 
specific  source- to- station  path.  One  of  the  surprising  results  of  this 
study,  as  well  as  related  ones  using  earthquake  data  for  South  America,  India 
and  the  Near  East,  is  that  the  Lg  coda  method  appears  to  give  reliable  average 
values  for  specific  source- to- station  paths.  This  observation  has  been  tested 
in  a  number  of  independent  ways,  such  as  using  Lg  amplitude  data  rather  than 
the  variation  of  predominant  wave  frequency  with  lapse  time  to  estimate  Q. 

Nevertheless,  it  is  instructive  to  look  at  the  consequences  of  errors  in 
Q  values.  If,  for  a  given  source  region,  the  Q  values  for  all  the  source-to- 
station  paths  are  overestimated,  the  resulting  m^Lg)  will  be  underestimated. 
Therefore,  the  absolute  m^Ig)  values  of  the  explosions  would  be  in  error, 
and  in  such  cases  these  m^(Lg)  values  could  not  be  used  to  estimate  bias 
(for  teleseismic  P  waves)  between  two  test  sites,  such  as  NTS  and  the  East 
Kazakh  sites.  However,  if  the  m^(Lg)  data  are  only  to  be  used  for  estimating 
yield  for  a  test  site  for  which  adequate  announced  yield  data  are  available 
for  calibration  purposes,  errors  in  rn^Lg)  would  not  affect  errors  in  yield 
estimates  for  that  site  as  long  as  the  yields  to  be  estimated  lie  in  the 
range  of  the  calibration  data.  Another  way  of  stating  this  is  to  note  that 
even  if  the  proposed  methodology  successfully  predicts  the  yields  of  NTS  ex¬ 
plosions  (those  whose  yields  are  classified),  there  is  no  assurance  that  the 
Q  values  employed  and  the  m^(Lg)  values  obtained  are  the  correct  values. 

After  a  number  of  iterations  our  best  Q  values  for  NTS  paths  to  DUG, 

BKS  and  TUC  ares  DUG,  Qq  =  155.  5  =  0.6j  BKS,  Qq  =  139,  T  =  0.6;  TUC,  Qq  = 

162,  3  =  0.6.  The  data  used  to  obtain  these  values  are  given  in  Tables  1  to  3, 


and  the  data  points  and  the  corresponding  curves  are  shown  in  Figure  1. 

Rondout  Associates  (1984)  also  applied  the  coda  Q  method  to  VWSSN 
statins  for  MTS  events.  They  obtained  smaller  1  values  and  larger  Qq  values 
than  those  given  by  Nuttli  (1983) »  namely t  DUG,  Qq  =  230,  5=  0.3;  BKS, 

Qq  =  225f  T =  0.2}  TOC,  Qq  =  300,  $  =  0.2.  As  a  result  their  estimates  of 
m^(lg)  for  NTS  events  axe  smaller  than  those  obtained  by 'Nuttli  (1983). 

Although  application  of  the  coda-Q  method  to  obtain  Qq  and  t  values  is 
fairly  straightforward,  there  are  two  complicating  factors  that  must  be 
watched  out  for.  First,  explosions  are  excellent  exciters  of  fundamental¬ 
mode  surface  waves  because  of  their  shallow  depth.  Depending  on  the  epi- 
central  distance,  the  fundamental-mode  waves  will  arrive  early  in  the  coda, 
about  one  to  three  minutes  after  the  onset  of  Lg.  They  will  have  larger  am¬ 
plitudes  and  longer  periods  than  the  scattered  coda  waves.  Therefore,  it  is 
easy  to  mistakenly  read  their  wave  frequency  rather  than  the  higher  frequency 
of  coda  waves  at  the  time  the  fundamental-mode  waves  are  arriving.  The  other 
problem  arises  near  the  end  of  the  coda.  In  this  study  it  has  been  observed 
that  beyond  some  particular  lapse  time,  whose  value  depends  on  the  size  of 
the  explosion,  the  majority  of  scattered  coda  waves  show  no  change  in  fre¬ 
quency  with  increasing  coda  lapse  time,  rather  than  continuing  to  decrease 
with  increasing  lapse  time.  One  has  to  carefully  search  that  part  of  the 
coda  for  the  relatively :few  and  somewhat  larger  amplitude  waves  that  give 
the  low  frequency  values.  If  not,  the  wave  frequencies  will  be  overestimated 
at  these  larger  lapse  times.  These  two  complicating  factors  can  lead  to 
smaller  estimates  of  S  and  consequently  larger  values  of  Qq. 

The  differences  between  the  Qq,T  values  for  paths  from  NTS  to  DUG,  BKS 
and  TUG  obtained  by  Nuttli  (1983)  and  Rondout  Associates  (1984)  lead  to  big 
differences  in  estimated  m^(Lg).  Table  4  illustrates  the  differences,  which 
amount  to  about  0.5  magnitude  units.  Comparison  with  m^(p)  values  published 
in  the  Bulletin  of  the  International  Seismological  Centre  lead  to  two  different 
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interpretations .  The  Rondout  Associates  m^(Lg)  values  are  about  0.2  of  a 
magnitude  unit  smaller  than  the  m^P)  values.  If  the  Rondaut  Associates 
m^(lg)  values  are  correct,  this  implies  that  Ig-wave  amplitudes  from  nuclear 
explosions  axe  smaller  than  from  earthquakes,  and  therefore  1-sec  period  Ig- 
wave  amplitudes  can  serve  as  a  good  discriminant  between  explosions  and 
earthquakes .  On  the  other  hand,  if  the  Nuttli  (1983)  m^(lg)  Va^-Ues  are 
correct,  ra^Ig)  for  NTS  is  on  average  about  0.3  magnitude  units  larger. than 
m^P-ISC).  If  these  m,o(Lg)  values  are  the  correct  ones,  it  follows  that  the 
teleseismic  P-wave  amplitudes  at  NTS  are  being  attenuated  by  at  least  0.3 
magnitude  units  in  the  asthenosphere  beneath  NTS.  (ihe  0.3  unit  value  would 
result  if  the  excitation  of  1-sec  period  Lg  amplitudes  for  earthquakes  and 
explosions  of  the  same  m^P)  values  is  the  same.  If  the  explosion  is  less 
efficient  in  generating  Lg  waves,  the  numerical  value  of  P-wave  attenuation  at 
NTS  would  be  even  larger  than  0.3  magnitude  units.) 

Earthquakes  near  NTS  offer  the  promise  of  deciding  between  the  above  two 
possibilities.  Dermengian  et  al  (1984)  gave  a  list  of  small  magnitude  earth¬ 
quakes  with  epicenters  near  NTS  that  occurred  between  1963  and  1968.  Of  the 
ten  earthquakes  in  their  list,  we  had  film  copies  of  three  in  our  seismogram 
library.  Table  5  presents  the  m^(P-NEIS)  and  the  m^(Lg)  values  obtained  by 
using  the  Nuttli  (1983)  estimates  of  and  T .  From  the  table  it  can  be 
seen  that  the  m^(Lg)  values  are  about  0.1  magnitude  units  larger  than  the 
m^(p)  values.  As  these  events  are  earthquakes,  m^(Lg)  is  expected  to  be 
larger  than  or  equal  to  mb(p),  i.e.  larger  if  there  is  anomalous  attenuation 
of  P-wave  amplitudes  in  the  asthenosphere  beneath  NTS,  and  equal  if  there  is 
no  such  anomalous  attenuation  of  P-wave  amplitudes.  On  average  for  the  three 
Nevada  earthquakes  m^Ig)  is  0.1  units  larger  than  mb(p)  if  the  Nuttli  (1983) 
Q,o  and  t  values  are  used.  If  the  Rondout  Associates  (1984)  and  T  values 
are  used,  the  m^Lg)  values  are  about  0.4  units  smaller  than  the  mb(p)  values 
for  the  earthquakes,  which  is  difficult  to  explain.  Although  it  would  be 
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desirable  to  have  more  data,  the  econcusion  Is'  that  the  Nuttli  (1983) 
values  of  and  X  are  more  applicable  to  the  NTS  region  than  the  Rondout 
Associates  (1984)  values. 


NTS  STUDIES 

Lg  amplitudes  were  measured  at  DUG,  BKS  and  TUC  for  all  NTS  events  of 
announced  yield  far  which  seismogram  copies  were  available.  The  events  were 
grouped  Into  two  simple  classifications,  alluvium  and  "hard  rock",  depending 
upon  the  characteristics  of  the  source  medium  as  described  in  Springer  and 
Kinnaman  (1971 »  1975)*  The  yields  also  were  taken  from  the  papers  of 
Springer  and  Kinnaman  (1971,  1975)*  P-wave  m^  values  were  obtained  from  the 
Bulletin  of  the  International  Sei sinological  Centre. 

The  m^(lg)  data  versus  log^Q  yield  are  plotted  in  Figure  2.  Over  the 
range  of  2  to  1000  kt  the  data  points  have  been  fitted  by  a  quadratic  curve, 
rather  than  a  straight  line.  If  the  Lg  spectrum  can  be  represented  by  a 
long-period  segment  of  zero  slope  and  a  short-period  segment  of  slope  two 
(the  period  at  which  the  two  segments  of  a  spectrum  intersect  is  called  the 
corner  period),  then  the  curve  of  Figure  2  should  consist  of  one  straight- 
line  segment  of  unit  slope  for  explosions  with  a  spectral  corner  period  of 
less  than  one  second  and  a  second  straight-line  segment  of  0.5  slope  for  the 
larger  explosions.  This  would  present  an  alternate  way  of  selecting  a 
mathematical  relation  to  fit  the  data.  The  equation  of  the  quadratic  curve, 
obtained  by  a  least-squares  fit  to  the  data,  is 

mb(lg)  =  3.868  +  1.181  log  Y  -  0.0938  (log  Y)2  (1) 

where  Y  is  the  announced  yield,  in  kilotons,  The  standard  deviation  of 
m^Lg)  is  0.060.  The  slope  of  the  curve  at  small  yields  (10  kt)  is  1.01, 
close  to  that  predicted  by  theory.  The  slope  at  the  largest  yields  shown  in 
the  figure  (1000  kt)  is  0.62. 

F*om  Figure  2  it  can  be  seen  that  the  two  granite  explosions  lie  above 


the  curve  (about  0.1  magnitude  unit)  and  the  dry  tuff  explosions  lie  slightly 
below  the  curve. 

The  alluvium  data  need  to  be  treated  separately,  as  they  clearly  form 
a  different  population  than  the  "hard  rock"  explosion  data.  Figure  3  shows 
the  data  and  a  curve  constructed  to  pass  through  the  data  and  to  have  the 
same  slopes  as  the  "hard  rock"  curve  at  yields  of  0.5  and  100  kt,  namely  I.05 
and  0.80,  respectively. 

In  addition  to  explosions  of  announced  yield,  m^(l«)  was  determined  also 
for  selected  NTS  events  whose  yield  values  remain  classified  information.  Of 
all  the  NTS  "hard  rook"  events  considered  in  this  study,  71  had  published 
mb(p-ISC)  values.  For  these  71  NTS  "hard  rock"  events  it  was  found  that 

m^Ig)  -  m^(P-lSC)  =  +0.31  -  0.02  (2) 

Assuming  that  m^(Lg)  is  equal  to  m^(P-ISG)  in  areas  where  there  is  no  anoma¬ 
lous  P-wave  amplitude  attenuation  in  the  asthenosphere  (e.g.  eastern  North 
America),  equation  (2)  tells  us  that  the  P-wave  amplitudes  are  reduced  by 
0.31  magnitude  units  in  the  asthenosphere  beneath  NTS.  If  Lg  amplitudes  are 
not  excited  as  strongly  for  explosions  as  for  earthquakes  of  the  same  mb(p) 
value,  the  amplitude  loss  of  P  waves  in  the  asthenosphere  beneath  NTS  would 
be  even  larger  than  0.31  magnitude  units. 

NON-NTS  EXPLOSIONS  IN  THE  UNITED  STATES 

There  are  a  few  explosions  of  announced  yield  in  the  United  States  not 
located  at  NTS.  Of  these,  SALMON  in  Mississippi,  GASBUGGY  in  New  Mexico 
and  RULISON  and  RIO  BLANCA  in  Colorado  were  selected  for  analysis.  Table  6 
gives  the  results,  where  the  NTS  hard-rock  yield  curve,  given  by  equation  (l), 
was  assumed  to  be  applicable  for  the  explosions  off  the  test  site.  In 
general,  the  agreement  between  announced  yield  and  the  n^Lg)  yield  estimate 
is  quite  good,  except  for  RIO  BLANCA,  an  event  of  three  30  kt  explosions, 
for  which  m^(lg)  underestimated  the  total  yield  by  27%.  However,  it  should 
be  noted  that  only  one  WWSSN  station  could  be  used  for  GASBUGGY,  RULISON  and 


and  RIO  BLANCA,  and  that  the  anelastic  attenuation  coefficient  value  from 
the  source  to  the  individual  station  was  obtained,  therefore,  only  from  the 
coda  of  one  seismogram.  Most  of  the  regional  stations  for  these  events  could 
not  be  used  because  the  lg  amplitudes  were  too  large,  with  the  trace  going 
off  the  seismogram. 

FRENCH  SAHARA  EXPLOSIONS 

Marshall  et  al  (1979)  gave  yields  of  several  French  Sahara  explosions. 
The  nearest  WWSSN  station  haveng  a  continental  path  was  HLW  in  Egypt,  at  a 
distance  beyond  3900  km.  Fortunately,  the  anelastic  attenuation  of  Lg  across 
Africa  is  low,  so  that  two  of  the  explosions,  SAPHTR  and  RUBIS,  had  measure- 
able  Ig-wave  amplitudes.  Figure  4  gives  an  example  of  Lg  waves  of  SAPHIR 
recorded  at  AAE  in  Ethiopa,  at  a  distance  of  394?  km.  The  onset  of  Lg  is 
sharp,  at  li*1  48®  55s*  The  maximum  sustained  amplitude  was  read  at 
llh  49®  08s.  The  seismogram  was  enlarged  to  approximately  six  times  its 
original  size,  and  the  Lg  waveform  was  traced  from  the  image  on  a  film 
viewer. 

Table  7  lists  the  data  for  RUBIS  and  SAPHIR.  Station  HLW  is  in  Egypt, 
and  SDB  in  Angola.  Both  shots  were  in  granite.  For  the  Nevada  explosions 
SHOAL  and  PILEDRIVER,  also  in  granite,  the  m^(lg)  was  about  0.15  units  larger 
than  for  the  typical  hard-rock  shot.  If  such  a  correction  were  applied  to 
the  Sahara  events,  the  estimated  yield  from  m^(lg)  for  RUBIS  would  be  46  kt 
and  for  SAPHIR  would  be  73  kt.  The  data  are  too  few  to  argue  whether  the 
French  Lg  data  are  better  satisfied  by  an  NTS  hard-rock  or  NTS  granite  m^(Lg) 
versus  yield  relation,  although  the  SAPHIR  data  would  argue  for  the  former, 
even  if  the  SDB  value  were  not  used.  The  conclusion  can  be  made,  however, 
that  for  the  limited  available  data,  the  NTS  hard-rock  m^(Lg)  versus  yield 
relation  will  give  French  Sahara  explosion  yields  no  more  than  50%  in  error. 


EXPLOSIONS  AT  EAST  KAZAKH  TEST  SITE 

Regional  Lg  waves  can  be  observed  at  VTWSSN  stations  in  southern  Asia 
and  Scandinavia  from  explosions  at  the  East  Kazakh  Test  Site  in  the  USSR. 
Table  8  lists  the  stations  utilized  in  the  study,  along  with  the  Qq  and  T 
values  obtained  by  applying  Herrmann's  (1980)  method  to  the  Lg  coda.  Also 
included  are  revised  Qq  values,  which  were  obtained  by  applying  station 
corrections  based  on  average  m^(lg)  values  for  the  events.  In  three  cases 
(KBL,  NIL  and  UNE)  the  revised  value  falls  outside  the  bounds  obtained 
by  application  of  the  coda  method.  In  all  three  cases  the  revised  Qq  is 
lower  than  the  lower-bound  value.  This  indicates  that  the  Lg  amplitudes 
at  these  stations  are  somewhat  smaller  than  expected. 

Figure  5  presents  m^(P-ISC)  versus  m^(Lg)  for  30  selected  Shagan  River 
events  at  the  East  Kazakh  Test  Site.  Contrary  to  NTS  events,  mb(P-ISC)  is 
larger  than  m^Lg)  for  Shagan  River  explosions.  The  average  difference  is 
0.10  units,  with  a  standard  deviation  of  0.14.  The  straight-line  curve 
drawn  in  the  figure  is  the  best-fitting  line  with  unit  slope.  Figure  6  shows 
similar  data  for  24  selected  Degelen  Mountain  explosions  at  the  East  Kazakh 
Test  Site.  Again  m^P-ISC)  is  larger  than  mb(Lg),  but  the  average  is  0.24 
units,  with  a  standard  deviation  of  0.14.  The  difference  in  average  values 
of  the  magnitude  difference  for  the  two  sites  at  East  Kazakh  is  surprising. 

It  is  difficult  to  explain,  if  it  is  real.  The  value  of  mb(P-ICC)  -  mb(Lg) 
for  PXLEDRIVER,  which  is  in  granite  at  NTS,  is  similar  to  the  values  for  the 
other  hard-rock  events  at  that  site.  Therefore,  the  different  source-rock 
types  at  Degelen  Mountain  and  Shagan  River  likely  are  not  the  explanation. 

The  difference  in  mb(P-ISC)  -  m^Lg)  values  at  the  two  East  Kazakh  sites, 
namely  0.24  -  0.10,  corresponds  to  about  a  38#  difference  in  estimated  yield 
values  based  on  mb(Lg). 


Tables  9  and  10  present  estimated  yields  of  selected  explosions  at  the 
Shagan  River  and  Degel en  Mountain  sites,  respectively.  They  are  not  intended 
to  be  taken  as  actual  estimates  of  the  true  yields  of  the  events,  but  rather 
to  demonstrate  the  methodology  for  using  m^(lg)  to  obtain  yield.  Before 
they  can  be  accepted,  it  must  be  demonstrated  that  the  NTS  hard-rock  relation, 
given  by  equation  (l),  applies  equally  well  to  Bast  Kazakh  explosions.  Also, 
the  effect,  if  any,  of  different  source  crustal  structures  on  the  excitation 
of  Lg  waves  needs  to  be  investigated.  Finally,  more  data  should  be  scruti¬ 
nized  to  obtain,  if  possible,  better  Qq  and  values  for  the  various  paths 
from  East  Kazakh  to  the  WVSSN  station. 

Marshall  et  al  (1979)  list  one  East  Kazakh  explosion  of  announced  yield. 
It  is  the  event  of  January  15,  1965  that  cratered  in  a  wet  sandstone,  for 
which  a  yield  of  125  kt  was  given.  The  estimated  yield  obtained  by  rn^Lg), 
namely  105  kt,  is  16%  less  than  the  announced  yield.  One  example  is  in¬ 
adequate  for  attempting  to  provide  explanations  for  the  difference.  It  even 
is  possible  that  the  relatively  close  agreement  is  fortuitous. 

The  data  of  Tables  9  and  equation  (2)  give  an  mb(p)  bias  between  NTS 
and  Shagan  River  of  0.41  magnitude  units,  subject  to  the  same  limitations 
as  the  yield  estimates  discussed  above.  This  bias  value  can  be  compared  with 
a  value  of  0.24  obtained  by  Murphy  and  Tzeng  (1982)  by  comparing  the  spectra 
of  Alaskan  earthquakes  as  obtained  from  seismograms  recorded  at  NTS  and  near 
the  Semipalatinsk  site,  and  of  0.25  obtained  by  Dermengian  et  al  (1984)  by 
means  of  m^/K  analysis  of  earthquakes  with  epicenters  near  the  Semipalatinsk 
site  and  NTS.  Sykes  and  Ciufentes  (1984)  used  20-sec  period  surface  waves 
from  NTS  events  to  obtain  an  Mg  versus  yield  calibration  curve,  which  they 
then  applied  to  East  Kazakh  explosions  to  estimate  yields.  They  also  calcu¬ 
lated  m^  values  for  the  NTS  and  East  Kazakh  events  which,  when  combined  with 
yield  values,  can  give  the  mfe  bias  between  the  two  test  sites.  Although 
they  did  not  explicitly  state  the  value  of  the  bias,  their  data  suggest  a 
value  of  at  least  0.35  magnitude  units. 


The  distribution  of  WWSSN  stations  is  fax  from  ideal  for  determining 
m^(lg)  and  for  estimating  yields  of  East  Kazakh  explosions.  The  southern 
Asian  stations,  which  are  closest  to  the  test  site,  are  located  in  or  beyond 
mountainous,  neotectonic  areas.  The  paths  to  the  Scandinavian  stations  are 
more  uniform,  but  the  distances  are  all  3500  km  and  greater.  Both  of  these 
problems  lead  to  scatter  in  rn^(Lg)  values.  A  0.1  unit  error  in  m^(Ig) 
translates  into  about  a  30#  error  in  estimated  yield.  The  ideal  data  would 
come  from  stations  in  the  interior  of  Asia,  at  distances  of  500  to  2000  km 
from  the  source. 
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TABLE  1 

WAVE  FREQUENCY,  f  ,  VERSUS  LAPSE  TIME,  t,  FOR  NTS  EVENTS  RECORDED  AT  BKS 
ON  WWSSIr  SHORT- PERIOD ,  VERTICAL-COMPONENT  SEISMOGRAMS 


Date  Event 

fp(Hz) 

t  (sec) 

02/08/6?  WARD 

4.77 

1.17 

171 

0.90 

200 

0.82 

213 

0.84 

260 

06/12/69  TAPPER 

4.85 

0.65 

0.98 

297 

213 

0.84 

261 

0.78 

110 

12/17/69  LCVAGE 

4.96 

0.73 

0.59 

351 

389 

0.56 

430 

08/31/67  DOOR  MIST 

5.05 

0.47 

1.07 

500 

188 

07/27/67  STANLEY 

5.19 

0.86 

1.01 

220 

199 

0.95 

214 

0.78 

233 

06/26/67  MIDI  MIST 

5.29 

0.62 

0.98 

290 

202 

0.84 

266 

01/30/69  VISE 

5.23 

O.65 

0.90 

320 

287 

0.76 

305 

0.49 

390 

05/10/67  MICKEY 

5.26 

0.42 

O.65 

440 

377 

O.59 

i>44 

0.50 

48? 

05/27/69  TORRIDO 

5.44 

0.43 

1.0? 

563 

185 

0.84 

205 

01/19/67  NASH 

5.51 

0.84 

0.90 

299 

240 

0.59 

355 

01/15/69  WINESKIN 

5-58 

0.53 

O.59 

395 

295 

0.49 

367 

0.43 

426 

0.37 

484 

10/29/69  CALABASH 

5.88 

0.29 

0.49 

498 

335 

0.42 

387 

09/27/67  ZAZA 

6.00 

0.39 

O.45 

437 

391 

0-39 

430 

0.33 

530 

..*  .*  ./  o  •_ 


.  w-. 
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TABLE  2 

.*  • 

WAVE  FREQUENCY,  fp,  VERSUS 

LAPSE  TIME,  t,  FOR  NTS 

EVENTS  RECORDED  AT 

DX 

3 

ON 

WWSSN  SHORT-PERIOD,  VERTICAL-COMPONENT  SEISMOGRAMS 

-  — ■* 

>*, 

Date 

Event 

mbC1*) 

fp(Hz) 

t(sec) 

>  *, 

»  *, 

02/23/67 

PERSIMMON 

4.2  5 

1.17 

165 

'•v> 

■> 

1.07 

179 

05/01/70 

0.98 

0.78 

224 

263 

■vv 

-  — 4 

BEEBALM 

4.73 

0.95 

235 

* 

0.59 

350 

. 

08/31/67 

DOOR  MIST 

5.05 

0.73 

320 

0.78 

357 

0.69 

388 

-  * 

2 

07/27/67 

STANLEY 

5-19 

0.73 

314 

j 

0.73 

369 

06/26/67 

0.59 

447 

MIDI  MIST 

5.20 

0.84 

302 

0.72 

373 

■ 

01/19/67 

0.51 

489 

». 

NASH 

5.51 

0.59 

374 

0.53 

433 

0.39 

496 

. 

0.45 

510 

02/23/67 

0.31 

574 

AGUE 

5-91 

0.98 

261 

a 

0.53 

455 

» 

0.39 

475 

l‘  m 

0.28 

870 

09/06/68 

NOGGIN 

6.08 

1.17 

150 

.*• 

0.98 

187 

•  ? . 

* 

0.90 

274 

s® 

0.69 

337 

0.45 

392 

06/06/73 

0.37 

4  65 

ALMENDRO 

6.31 

0.37 

500 

0.34 

533 

0.28 

574 

i 

0.24 

593 

09/16/69 

JORUM 

6.52 

0.37 

489 

v' 

0.32 

508 

SV 

0.24 

592 

0.26 

610 

*  * 

0.22 

680 

—  4 

►  m'» 

*  .  % 

i  _ 

»  /  • 

,*‘  *  < 

•!v 

\  .\i 

*  » 

7 

'  i 

.*•  - 

>  * 

.  -  -  V.  ‘  .  .  ■  .  •  » 

i  . 

:.-..\v.\-.v  v. 
_v. 
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TABLE  3 

WAVE  FREQUENCY,  fp,  VERSUS  LAPSE  TIME,  t,  FOR  NTS  EVENTS  RECORDED  AT  TUC 
ON  WWSS1T  SHORT-PERIOD,  VERTICAL-COMPONENT  SEISMOGRAMS 


Date 

Event 

“bt1*) 

fp(Hz) 

t(sec) 

02/23/67 

PERSIMMON 

4.25 

0.99 

233 

0.90 

296 

02/08/6? 

WARD 

4.77 

0.84 

0.84 

316 

285 

0.78 

310 

06/29/67 

UMBER 

4.82 

0.73 

0.82 

355 

316 

O.65 

410 

08/18/67 

BORDEAUX 

5.01 

0.59 

0.73 

438 

274 

0.62 

382 

07/27/67 

STANLEY 

5.19 

O.51 

0,84 

4  55 
297 

O.65 

372 

0.49 

491 

06/26/6 7 

MIDI  MIST 

5.20 

0.42 

0.67 

555 

390 

0.49 

425 

01/30/69 

VISE 

5-2 3 

0.43 

0.53 

450 

393 

0.45 

426 

05/10/67 

MICKEY 

5-26 

0.45 

O.65 

475 

377 

0.59 

444 

0.50 

487 

01/19/67 

NASH 

5.51 

0.43 

0-53 

563 

440 

0.51 

487 

05/07/69 

PURSE 

5-98 

0.36 

0.25 

573 

622 

0.23 

668 

10/08/69 

PIPKIN 

6.06 

0.20 

0.53 

793 

427 

0.42 

433 

0.28 

510 

0.39 

535 

0.33 

625 

0.28 

665 

09/16/69 

JORUM 

6.52 

0.23 

0.21 

753 

701 

0.20 

1080 

0.20 

1140 

01/03/76 

MUENSTER 

6.53 

0.20 

0.20 

1290 

708 

0.17 

770 

0.19 

890 

TABLE  5 


iab(P-NEIS)  AND  m^(Lg)  VALUES  FOR  EARTHQUAKES  NEAR  NTS 


Date 

08/21/64 

Origin  Time 

22-03-51.65 

Latitude  (°N) 
37-0 

Longitude  (°W) 
115.1 

Vp) 

3.8 

11/17/65 

09-41-28; 3 

37.6 

115.2 

3.7 

04/06/66 

17-56-32.1 

37.28 

115-38 

4.1 

4.1 
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TABI£  5 


ESTIMATED  m.(lg)  FOR  SEIECTED  NTS  EVENTS  USING  DIFFERENT 
ATTENUATION  VALUES 


Date 

Event 

^(P-TSC) 

Rondout  expat' 1 

Ronaout  linear 

nk/Lg) 

Nuttli 

09/13/63 

BILBY 

— 

5.59 

5.5? 

6.1 

02/24/66 

REX 

5.0 

4.86 

4.78 

5.28 

05/06/66 

CHARTREUSE 

5.4 

5.24 

5. 22 

5-76 

05/20/67 

COMMODORE 

5.8 

5.52 

5.50 

6.04 

05/23/67 

SCOTCH 

3.7 

5.44 

5.44 

6.00 

12/19/68 

BENHAM 

6.3 

6.13 

6.01 

6.65 

04/26/73 

STARVQRT 

5.6 

5-35 

5-33 

5.85 

TABUS  6 

ESTIMATES  OF  YEILD  OF  NON-NTS  EXPLOSIONS  IN  CONTINENTAL  U.  S. 


_  .  _  .  Stations 

Event  Date  Used  ^(i*)  a 

SALMON  10/22/64  BLA,  DAL,  4.66 

FLO,  OXF,  ROD 

GASBUGGY  12/10/67  DUG  5.38 

RULISON  09/10/69  TUC  5.56 

RIO  BLANCA  O5/17/73  TUC  5.71 

♦Uses  NTS  hard-rock  relation,  equation  (l). 


mb(Lg)  m^(P-ISC) 


Announced 
yield  (kt) 


12/10/67 

09/10/69 


Yield  from 
mb(Lg)  (kt) 


l,  DAL,  4.66 

4.6 

5.3 

5 

OXF,  RCD 

DUG  5*38 

4.8 

29 

28 

TUC  5.56 

5-0 

40 

45 

TUC  5.71 

5.1 

90 

66 

TABLE  7 

FRENCH  SAHARA  EXPLOSIONS 


Event 

Date 

Station 

L  (km) 

Qo 

Y 

®b(Dg) 

Yield 

mbCDg 

RUBIS 

1 

CM 

H 

AAE 

3947 

750 

0.3 

5.72 

52  ’ 

68 

SAPHTR 

02/27/65 

AAE 

3947 

750 

0.3 

3-96 

(135) 

HLW 

3146 

700 

0.3 

5-92 

(120) 

SDB 

4414 

720 

0.4 

5.79 

(82) 

average 

5.89 

120 

111 

♦Uses  NTS  hard-rock  relation,  equation  (l). 


TABLE  8 

ATTENUATION  PARAMETERS  FOR  PATHS  FROM  EAST  KAZAKH 
TEST  SITE  TO  SELECTED  WWSSN  STATIONS 


Station 

Approximate 
Distance  (km) 

Q 

Coda  Method 

Upper  &  Lower  Bounds 

Y 

Revised 

Qo 

KBL 

1850 

360 

(340,  380) 

0.6 

336 

KEV 

3500 

580 

(550,  620) 

0.4 

554 

KON 

4350 

700 

(650,  760) 

0.4 

700 

MHI 

2100 

380 

(360,  410) 

0.5 

380 

NDI 

2350 

300 

(270,  320) 

0.6 

312 

NIL 

1850 

380 

(370,  400) 

0.6 

354 

NUR 

3500 

580 

(520,  650) 

0.4 

580 

QUE 

2150 

300 

(280,  330) 

0.6 

318 

SHL 

2950 

340 

(325,  350) 

0.6 

340 

UME 

3700 

620 

(610,  64o) 

0.4 

591 

TABUS  9 


TENTATIVE  YIELD  ESTIMATES  FOR  SELECTED  SHAGAN  RIVER 
EXPLOSIONS  OBTAINED  BBOM  m^Lg)  VALUES 

Notes  The  estimated,  yields  axe  given  only  for  a  demonstration  of  the  methodology, 
and  should  not  be  interpreted  or  used  as  measures  of  actual  yields.  A 
number  of  presently  unverified  assumptions  were  made  in  obtaining  the 
estimated  yields. 


Date 

01/15/65 

11/02/72 

12/10/72 

07/23/73 

05/31/74 

10/16/74 

12/27/74 

04/27/75 

10/29/75 

12/25/75 

04/21/76 

06/09/76 

07/04/76 

08/28/76 

12/07/76 

05/29/77 

06/29/7? 

09/05/77 


06/11/78 

07/05/78 

08/29/78 

09/15/78 

11/04/78 

04/05/79 

06/23/79 

07/07/79 

08/04/79 

08/18/79 

09/14/80 


5.8 

6.1 

6.0 

6.1 

5-9 

5-.5 

5.6 

5.6 

5.8 
5-7 
5.3 
50 
5-8 
5-8 

5.9 
5-8 
5-3 
5-8 
6.0 
5.9 

5.8 

5.9 
6.0 
5-6 
5.8 

6.2 

5.8 

6.1 

6.1 

6.2 


Number  of 
Stations 

Estimated., 
Yield  (kt) 

5.87 

2 

105 

6.04 

1 

170 

6.09 

2 

200 

6.13 

6 

229 

5.68 

4 

62 

5.26 

3 

20 

5.69 

6 

63 

5-47 

4 

35 

5.45 

2 

33 

5.83 

4 

93 

5.19 

1 

17 

5.27 

3 

21 

5-90 

2 

112 

5.60 

5 

49 

5. 71 

2 

66 

5-5 8 

5 

47 

5-15 

3 

16 

5.51 

2 

39 

5-71 

5 

66 

5-75 

4 

74 

5.67 

1 

59 

5.80 

3 

87 

5- 87 

3 

105 

5-57 

2 

46 

5.84 

2 

96 

5.92 

2 

120 

5-87 

3 

105 

6.01 

4 

157 

6.03 

5 

166 

6.09 

5 

200 

Announced-, 
Yield  (kt) 

125 


♦m^CP)  values  from  International  Seisraological  Centre  Bulletin. 
**  Uses  NTS  hard-rock  relation,  equation  (l). 

♦♦♦Announced  yield  value  given  in  Marshall  et  al  (1979). 
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TABLE  10 

TENTATIVE  YIELD  ESTIMATES  FOR  SELECTED  DEGEIEN 
MOUNTAIN  EXPLOSIONS  OBTAINED  FROM  ^(Lg)  VALUES 

Note*  The  estimated  yields  are  given  only  for  a  demonstration  of  the  methodology, 
and  should  not  be  interpreted  or  used  as  measures  of  actual  yields.  A 
number  of  presently  unverified  assumptions  were  made  in  obtaining  the 
estimated  yields. 


Date 

“b(* **># 

Number  of 
Stations 

Estimated 
Yield  (kt) 

04/19/73 

5.4 

5.12 

2 

15 

07/10/73 

5.2 

5.14 

2 

15 

10/26/73 

5-2 

4.88 

1 

8.3 

01/30/74 

5.4 

5.11 

2 

14 

05/16/74 

5.2 

4.89 

2 

8.2 

07/10/74 

5-2 

4.75 

2 

6.2 

09/13/74 

5-2 

4.74 

1 

6.2 

12/16/74 

4.8 

^.15 

2 

1.7 

02/20/75 

5-7 

5.40 

3 

30 

03/11/75 

5.4 

5.23 

2 

19 

06/08/75 

5-5 

5.04 

2 

12 

08/07/75 

5-2 

4.90 

2 

8.7 

01/15/76 

5-2 

4.92 

2 

9.3 

04/21/76 

5-1 

5-21 

1 

18 

05/19/76 

5.0 

4.92 

2 

9.3 

07/23/76 

5.1 

4.96 

2 

10 

12/30/76 

5.2 

4.99 

1 

11 

03/29/77 

5-4 

5.36 

4 

26 

04/25/77 

5*1 

4.94 

2 

9-5 

07/30/77 

5-1 

4.76 

3 

6.3 

12/26/77 

4.9 

4.73 

2 

6.0 

03/19/78 

5-2 

4.90 

2 

8.7 

03/26/78 

5.6 

5.39 

5 

29 

04/22/78 

5-3 

5-19 

4 

17 

11/29/78 

6.0 

5-74 

4 

72 

*m^(P)  values  from  International  Seismological  Centre  Bulletin. 

**  Uses  NTS  hard- rock  relation,  equation  (l). 
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mb(P) 

m^P)  Versus  mb(Lg)  Values  for  Explosions  at  Degelen  Mountain 
East  Kazakh. 


CODA-Q  STUDIES  FOR  THE  INDIAN  SUBCONTINENT 
Viera  M.  John  and  Otto  V.  Nuttli 

The  coda-Q  method  as  developed  by  Aki  and  Chouet  (1975)  and  extended 
by  Herrmann  (1980)  was  applied  to  the  Lg  coda  recorded  by  short-period, 
vertical-component  WVSSN  seismograms  for  the  Indian  stations  KOD,  NDI, 

POO,  SHL  and  the  Pakistani  stations  NIL,  QUE.  All  six  stations  are 
located  in  the  Indian  subcontinent. 

Figure  1  shows  the  location  of  the  stations  and  the  earthquake  epi¬ 
centers  used  in  the  study.  The  epicenters  near  POO  are  the  Koyna  earthquake 
of  1967  and  its  aftershocks.  A  total  of  134  earthquakes  was  studied,  of 
which  7  were  near  Koyna  Dam. 

Figure  2  shows  the  results  of  the  study.  The  first  number  in  paren¬ 
theses  is  and  the  second  is  f  ,  where 

Q(f)  =  Q0  f  * 

and  f  is  frequency  and  Qq  is  the  value  of  Q  at  a  frequency  of  1  Hz. 

From  Figure  2  it  is  apparent  that  the  southern  shield  and  upland  area 
of  India  is  marked  by  high  Qq  and  relatively  low  values.  The  Qq  values 
are  three  or  more  times  larger  than  in  the  mountainous  areas  to  the  north. 
The  numbers  are  placed  on  the  map  about  midway  between  the  earthquake  and 
the  station. 

Figure  3  is  a  generalized  geological  and  structural  map  of  the  area. 
Although  the  shield  area  exhibits  structural  features,  it  is  of  relatively 
low  seismicity,  similar  to  eastern  North  America,  eastern  South  America  and 
central  and  northern  Asia. 

Figures  4  to  8  show  some  typical  examples  of  the  data  obtained.  Figure 
4  shows  the  plot  of  KOD  data,  for  the  earthquakes  which  are  to  the  northeast 
of  station  SHL.  Figure  5  presents  the  NIL  data  as  obtained  from  earthquakes 
with  epicenters  near  39°  N,  73°  E.  Figure  6  shows  the  data  at  POO  for  the 
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nearby  196?  Koyna  earthquakes.  Figure  7  presents  QUE  data  for  earthquakes 
near  39°  N,  73°  E.  The  Qq  value  of  265  is  noticeably  smaller  than  the  value 
of  640  obtained  from  QUE  seismograms  for  the  Koyna  epicenters  at  17°  N,  73°  E, 
The  phenomenon  exhibited  in  Figures  7  and  8,  namely  different  Qq  and  I  values 
for  earthquakes  from  different  source  regions  recorded  at  the  same  station, 
also  was  observed  at  NDI,  NIL,  POO  and  SHL. 

As  a  check  on  the  coda  Q  values,  Lg  amplitudes  were  observed  at  the 
same  station  for  a  number  of  earthquakes  at  different  distances  and  equalized 
to  m^  =  5*0  values.  The  resulting  Lg  amplitude  versus  distance  data  were 
plotted  and  fitted  by  a  curve  in  which  f ,  the  coefficient  of  anelastic 
attenuation,  is  a  parameter.  The  quantity  If  is  related  to  Q  by 

*(f)  =  rrf/u(f)Q(f) 

where  U  is  group  velocity.  The  value  of  X  obtained  by  amplitude  fall-off 
for  waves  of  a  given  frequency  was  compared  to  that  obtained  from  .the  Qq  and 
%  values  by  means  of  the  coda  analysis.  As  an  example,  for  station  NIL  the 
Lg  amplitude  attenuation  method  gave  a  Qq  of  300,  compared  to  a  value  of  330 
obtained  by  the  coda-Q  method.  The  agreement  is  considered  sufficiently 
good  to  justify  the  previously  arrived  at  conclusion  (Herrmann,  1980)  that 
coda-Q  methods  give  Qq  values  which  are  equal  to  the  apparent  Qq  values  of 
Lg  waves. 
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Figure  5*  Plot  of  fp(Hz)  versus  t  (sec)  for  NIL-  Theoretical  Curve  for 
Qq  =  2q5>  1  =  0*3*  The  Epicenter  Location  of  the  Earthquakes 
Plotted  is  39°  N,  73°  E. 
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Figure  8.  Plot  of  f  (Hz)  versus  t  (sec)  for  QUE.  Theoretical  Curve  for 


Qq  =  640,  ?  =  0.3.  The  Epicenter  Location  of  the  Earthquakes 
Plotted  is  17°  N,  73°  E. 


Attenuation  of  High  Frequency  Earthquake 


waves  in  South  America 
Mehdi  M.  Raoof  and  Otto  W.  Nuttli 
ABSTRACT 

Earlier  attenuation  studies  for  the  South  American  con¬ 
tinent  indicate  that  there  is  low  attenuation  for  Sn  and  Lg  waves 
in  the  shield  region  east  of  the  Andes,  whereas  in  the  west  of 
South  America,  with  some  exceptions,  there  is  high  attenuation 
for  Sn  and  Lg  waves. 

However,  these  studies  were  non-quantitative.  In  this  study, 

(1-Hz  values)  for  Lg  waves  for  South  America  are  presented, 

based  on  a  scattering  model  of  Aki  (1969)  as  extended  by 

Herrmann  (1980)  for  the  coda  waves  of  shallow  local  and  near- 

regional  earthquakes.  The  results  of  the  coda  Q  method  are 

compared  with  those  obtained  using  Nuttli’s  (1973)  method. 

These  coda-Q,  values  are  in  good  agreement  with  the  apparent  Qe 

of  Lg  waves  obtained  by  the  latter  method.  The  data  were 

obtained  from  over  100  local  and  regional  earthquakes  recorded 

by  12  WWSSN  stations  throughout  continental  South  America. 

These  data  provided  a  range  of  frequencies  from  0.4  to  1.4  Hz. 

Frequency  dependence  of  Q  was  investigated  for  the  observed 

range  of  frequencies  by  assuming  Q  =  Q0  (///„)*.  The  observed 

data  indicate  that  the  tectonic  region  of  western  South  America 

Department  of  Earth  and  Atmospheric  Sciences,  Saint  Louis  University,  St.  Louis, 
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is  characterized  by  low  Qa  and  a  large  value  of  the  frequency 
dependent  factor,  f,  with  values  ranging  from  150  to  350  and  0.4 
to  0.7,  respectively.  values  increase  in  the  shield  region  east 
of  the  Andes  but  frequency  dependence  decreases.  Average  cru¬ 
stal  Q,  values  obtained  for  north  and  central  Argentina  range 
from  420  to  580  and  f  ranges  between  0.2  and  0.3.  The  Q0  values 
are  larger  in  the  Brazil  region,  ranging  from  580  to  980,  with  f 
varying  from  0.0  to  0.2.  In  the  lower  attenuation  region  of 
eastern  South  America,  higher  values  of  attenuation  correlate 
with  greater  thickness  of  the  sedimentary  layers. 

INTRODUCTION 

The  study  of  the  anelasticity  of  various  regions  of  the  earth 
is  of  major  importance  for  earthquake  seismology.  In  the  con¬ 
tinental  United  States  and  in  Eurasia,  lateral  variations  in 
attenuation  and  frequency  dependence  of  the  quality  factor  Q 
have  been  observed.  In  particular,  the  short-period  surface 
waves  used  to  study  the  earth’s  crust  suffer  strongly  from  lateral 
heterogeneities  in  the  crust.  Therefore,  a  practical  analysis 
method  is  required  to  deal  with  the  problem  of  a  laterally 
heterogeneous  earth. 

To  deal  with  this  problem,  Aki  (1969)  proposed  a  statistical 
treatment  for  the  waves  which  are  backscattered  due  to  lateral 
heterogeneity.  These  waves  are  found  in  the  coda  of  local  earth¬ 
quakes.  He  suggested  that  the  seismic  coda  waves  of  local  earth¬ 
quakes  are  backscattered  waves  from  numerous  randomly  distri- 


buted  heterogeneities  in  the  earth.  Using  coda  waves  of  local 
earthquakes,  Aki  and  Chouet  (1975),  Aki  (1980a,  b),  Herrmann 
(1980),  Rautian  and  Khalturin  (1978),  and  Singh  (1961)  investi¬ 
gated  the  quality  factor  (Q)  of  seismic  waves  for  frequencies 
greater  than  0.1  Hz.  Herrmann’s  (1980)  Q  values  for  selected 
United  States  locations  are  consistent  with  Q  values  calculated 
from  the  anelastic  attenuation  coefficient  (7)  by  Nuttli  (1973). 

Because  coda  Q  values  are  consistent  with  apparent  Q  values  of 
Lg  waves  in  the  crust,  they  give  some  measure,  at  least  in  a  rela¬ 
tive  sense,  of  the  Q  structure  of  the  crust. 

The  purpose  of  this  study  is  to  obtain  apparent  Q  values  for 
high  frequency  Lg  waves  (approximately  1  Hz)  traveling  in  the 
South  American  crust  by  using  Herrmann's  (1980)  coda  Q 
method.  These  results  will  be  compared  with  those  obtained  by 
using  a  different  method,  the  amplitude  decrease  with  distance, 
which  was  used  by  Nuttli  (1973)  for  eastern  North  America. 

Previous  anelastic  attenuation  studies  have  been  primarily 
concerned  with  the  mantle  beneath  the  South  American  con¬ 
tinent.  Body  wave  attenuation  studies  by  Sacks  (1971).  Is  inks 
and  Barazangi  (1980),  Chinn  et  a l.  (1980)  and  others  indicate 
that  there  are  some  high  attenuation  zones  in  the  west  of  S  jut.h 
America. 

Crustal  studies  in  South  America  indicate  great  change  in 
the  crustal  thickness  from  the  coast  towards  the  central  \rdcs 
As  can  be  seen  in  Figure  1.  the  thickness  of  the  crust  changes 
from  about  7  km  sea-ward  of  the  trench  in  western  South  Amer- 
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iea  (Whitsett,  1975)  to  a  maximum  of  70  km  under  the  western 
high  plateau  (Altiplano)  of  the  Andes  (James,  1971).  The  crust 
then  thins  toward  the  interior  of  the  continent,  with  a  value  of 
25-40  km  to  the  east  of  the  eastern  cordillera.  Such  variations  in 
the  thickness  have  a  great  influence  on  the  propagation  of 
phases  such  as  I4.  Since  there  has  been  a  limited  number  of  Lg 
attenuation  studies  for  the  South  American  continent,  with  no 
Qi§  values  reported,  this  study  was  undertaken  to  provide  a  gen¬ 
eralized  Q0  map  for  Lg  waves  of  about  1-Hz  frequency  for  the 
South  American  continent,  as  well  as  a  contoured  Qa  map  of 
South  America. 


THEORY  AND  METHOD 

Aki  (1969)  observed  that  short-period  surface  waves,  which 
were  used  in  crustal  studies,  suffered  strongly  from  later  il 
heterogeneities  in  the  crust.  Therefore,  he  proposed  a  statistical 
treatment  to  deal  with  the  heterogeneity  of  the  earth.  »rd 
defined  the  seismic  coda  waves  of  local  earthquakes  as  backs  ut¬ 
tered  waves  from  numerous  randomly  distributed  het.-rj- 
geneities  in  the  earth.  Such  waves  can  be  found  even  near  K  e 
epicenter  after  the  passage  of  all  the  primary  waves  Therefore, 
the  later  portions  of  a  seismogram  should  be  the  rc>uli  if 
averaging  over  many  samples  of  heterogeneities 

Some  important  characteristics  of  the  coda  are  its  indepen¬ 
dence  of  epicentral  distance  and  azimuth  and  the  fact  that  -is  :- 
ful  properties  of  it  can  be  simply  measured  from  scismogr  in  s 


without  any  sophisticated  calculations. 

For  the  average  peak-to-peak  amplitude  observed  around 
time  t,  Aki  (1969)  proposed  the  following  relation: 

JiL.  (1) 

where  I(/p)  is  the  instrument  magnification  at  predominant  fre¬ 
quency  observed  at  time  t,  Q  is  the  anelastic  attenuation  quality 
factory  for  the  surface  wave  comprising  the  coda,  Afe  is  seismic 
moment  andB(fp)  is  defined  as 

(2) 

where  N(rB)  is  the  number  of  scatterers  within  a  radius  of  r0. 
9»b(/,,  rB)  is  a  function  of  frequency  only,  defined  by 

*,(vj\r)s\F(w\r)/Afaexp(.wt/2Q)  (3) 

where  F(w|r)  expresses  the  expected  absolute  value  of  the 
Fourier  transform  of  secondary  waves  coming  from  a  scattorer 
at  distance  r.  So  pB(wjr)  expresses  the  excitation  of  the  secon¬ 
dary  scattered  waves.  The  unit  of  pB(w|r)  is  cm  per  Hz  per  dyne 
cm. 

Equation  (1)  was  developed  by  Aki  to  analyze  earthquakes  n 
central  California  for  a  particular  value  of  Q.  Herrmann  (1183) 
expressed  equation  (l)  in  terms  of  another  time  variable  ‘  *  = 
t/Q,  to  be  applicable  in  other  areas,  as  follows: 


Q~xy2Sf,B  (fp  )C  (fp.t') 


U) 
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whore  the  term  C (/, ,  t*)  is  called  the  coda  shape  function  and  is 
equal  to 


*f. 


dt * 


1/4 


(5) 


For  the  application  of  equation  (5),  there  should  be  a  change 
of  predominant  frequency  of  the  coda  with  time,  dfp/ dt*. 

In  Table  1,  the  predominant  frequency  fp  and  the  coda 
•hape  function  C(/p,  t*)  are  given  as  a  function  of  t*  for  the 
short  period  WWSSN  instrument  system  with  a  peak 
magnification  at  1.35  Hz,  and  a  relative  magnification  of  1.0  at  1 
Hz.  Since  t*  *  t/Q,  if  predominant  frequency  fp  and  coda  shape 
are  plotted  as  a  function  of  t*.  they  result  in  theoretical  master 
curves  for  determination  of  Q  for  different  regions. 

By  assuming  a  source-spectrum  corner  frequency  greater 
than  the  peak  frequency  of  the  instrument  response,  the  effect 
of  source  spectrum  can  be  ignored  in  the  above  calculation.  Fre¬ 
quency  dependence  of  Q  can  be  investigated  by  assuming  that 


9(/)  =  9.  (///.)<  (6) 

assigning  a  value  to  ?,  and  then  plotting  fp  as  a  function  of  t* 
Master  curves  for  different  values  of  f  are  given  in  Figure  2. 

Methods  of  Q  Estimation 

There  are  several  methods  that  have  been  proposed  for 
estimating  Q.  some  of  which  are:  /,  versus  t  and  coda  shape 


methods  (Q  coda)  by  Herrmann  (1980),  and  Q  inferences  from 
gamma  )  by  Nuttli  (1973). 

In  this  study  the  f9  versus  t  method  will  be  used  first  by 
plotting  the  predominant  frequency  of  the  coda  as  a  function  of 
time  after  the  origin,  including  the  Lg  arrival.  The  predominant 
frequency  can  be  calculated  by  simply  counting  the  number  of 
zero  crossings  of  the  seismic  trace  within  a  given  time  interval. 
We  use  a  smaller  window  or  time  interval  shortly  after  the  Lg 
arrival  but  a  larger  one  for  later  parts  of  the  coda  because  of  the 
decreasing  frequency  content  in  the  coda.  The  number  of  zero 
crossings  is  divided  by  twice  the  window  length  to  obtain  fp. 
Next,  the  plot  of  observational  data  is  superimposed  on  the  mas¬ 
ter  curves  of  Figure  2- 

We  obtain  the  frequency  dependence  parameter  f  from  the 
best>fit  curve,  and  Qp  is  obtained  from  that  value  of  t  which  maps 
into  t*  *  1,  since  t*  »  t/Q. 

Figure  3  is  the  plot  of  predominant  frequency  fp  (Hz)  as  i 
function  of  time  t  (sec)  after  the  origin  at  the  Peldehue  (P£L) 
station.  The  theoretical  curve  for  Q,  equal  to  260  and  <■  =  0  •;  is 
also  shown. 

Since  the  /,  versus  t  data  plot  shows  considerable  scatter- 

m 

mg  for  most  of  the  stations  in  South  America,  which  sometir  .es 
makes  it  impossible  to  estimate  the  frequency  dependence  fac¬ 
tor  and  thus  the  best  value  of  Q,,  a  test  of  the  validity  of  1  he 
results  obtained  from  the  /,  versus  t  method  will  be  made  by 
using  also  Nuttli's  (1973,  1980)  method. 


The  amplitude  of  dispersed  surface  waves  in  the  time 
domain  is  given  by  Ewing  et  aL  (1957)  as 

A  \=A  tlA~l/i(R  a3mAa)“l/2exp(-7A)  (7) 

and  for  the  Airy  phases,  which  Nuttli  (1973)  concluded  gave  the 
best  fit  to  Lg  amplitude  attenuation,  as 

Aa=i40jA"l/,(j[?03inAa)"1/Jerp(-7A)  (8) 

where  A i  and  A2  ere  observed  amplitudes  at  an  epicentral  dis¬ 
tance  A,  Aqi  and  A  &  are  constants,  R0  is  the  radius  of  the  earth 
model,  and  7  is  the  coefficient  of  anelastic  attenuation.  Nuttli 
(1973)  used  the  fact  that  7  is  related  to  Q,  frequency  f,  and  group 
velocity,  U,  by 


tJSL 

QU 


O) 


The  group  velocity.  U,  for  the  Lg  phase  is  assumed  to  be  3.5 
km/sec.  The  term  A""  represents  a  decrease  in  amplitude  due 
to  dispersion  effects,  with  n  equal  to  1/2  for  ordinary  disper  ;ed 
surface  waves  and  1/3  for  an  Airy  phase.  The  next  term 
(i?,su»A,)”,/i  represents  the  amplitude  decrease  due  to  geometri¬ 
cal  spreading,  where  A  is  the  epicentral  distance  in  degrees  .md 
R9  is  the  radius  of  the  earth.  The  last  term  in  both  equations, 
exp  (-7A),  represents  attenuation  due  to  absorption.  Equation 
(9)  can  be  written  as 


For  a  given  amplitude  A  at  an  epicentral  distance  A  and  an  ampli¬ 


tude  A  jo  at  a  distance  of  10  km,  we  have 


A,0(/)=A(A./) 


L(fcm?  »[  «n(A»)  1  a 

P  10  «a<-l°-) 

l  '  111.1 '] 


»[y(A— 10)] 


where  f  is  Lg  wave  frequency. 


To  study  the  attenuation  of  Lg  amplitude  with  distance, 
Nuttli  (1973)  used  equation  (10)  to  get  a  set  of  master  curves, 
consisting  of  log-log  plots  of  amplitude  (A)  versus  distance  (A), 
where  each  curve  represents  a  different  value  of  y  (Figure  4).  If 
the  Lg  amplitude  is  plotted  as  a  function  of  epicentral  distance 
for  an  event  on  the  same  log-log  scale  as  the  master  curves,  we 
can  determine  the  value  of  y  by  fitting  the  curves  to  the  data, 
such  as  was  done  by  Nuttli  (1973),  Nuttli  and  Dwyer  (1978).  and 
Ballinger  (1979).  If  y  is  determined  in  this  manner  for  1-Hz  Lg 
waves,  Q,  can  be  obtained  from  equation  (9).  By  applying  the 
same  procedure  to  Lg  waves  of  various  frequencies,  f  also  can  be 
obtained  from  the  variation  of  y  and  Q  with  frequency. 

After  7  (1  Hz)  has  been  determined,  equation  (11)  can  be 
used  to  extrapolate  A(A)  to  obtain  A  jo-  For  eastern  North  Amer¬ 
ica,  Nuttli  (1973)  found  that  4l0  (1  Hz)  =110  microns  for  m*  = 
5.0,  where  m*  is  the  teleseismic  P-wave  magnitude  Making  this 
assumption 


log  A  JO  =  log  Am  +  (m»  +•  Am*  -  5.0) 
where  AM  =  A  to  for  m*  =  5.0.  Am* 
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is  the  m*  correction  parameter.  Therefore,  for  a  given  earth¬ 
quake,  if  m*  is  known  from  teleseismic  P-wave  amplitudes,  Ai0 
can  be  calculated.  If  the  excitation  of  Lg  waves  is  independent  of 
crustal  structure,  as  shall  be  assumed,  AM  ~  110  microns  every¬ 
where.  However,  for  the  same  strength  of  source  excitation,  m6 
varies  geographically  because  of  variation  in  the  anelastic 
attenuation  of  P-wave  amplitudes  in  the  asthenosphere.  This  is 
accounted  for  by  adding  aimt  term  to  the  value  of  m*  in  equa¬ 
tion  (12).  If  the  Am*  value  can  be  found  for  a  geographic  region, 
equation  (12)  will  give  a  1-Hz,  vertical  component  Lg  amplitude 
at  10  km  distance.  This  value  can  be  compared  with  that 
obtained  from  equation  (11),  and  the  best  set  of  Qa.  f  values  (to 
make  the  two  A «  values  agree)  of  the  permissible  set  of  values 
from  the  coda  Q  method  can  be  chosen.  This  procedure  helps  to 
reduce  the  uncertainty  in  estimating  f  when  the  range  of  fp.  t 
values  is  not  large  enough,  or  when  there  is  scatter  in  the  data 
paints,  so  that  Q,  and  f  values  cannot  be  determined  uniquely 
from  the  coda  data. 


DATA  ANALYSIS 

In  this  study  two  main  difficulties  were  encountered  They 
are  the  scarcity  of  data  due  to  low  instrument  magnification  and 
the  scatter  of  the  data.  Therefore  several  thousands  of  shallow 
and  deep  events  recorded  by  12  WTfSSN  stations  were  tested, 
from  which  over  240  local  and  regional  earthquakes,  occurring 
mostly  in  the  west  of  the  South  American  continent,  were 


.vVvl' 
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selected  for  the  determination  of  Q.  It  was  found  that  the  deep 
earthquakes  didn’t  give  results  consistent  with  shallow  earth* 
quakes.  Therefore,  the  deep  earthquakes  were  removed  from 
the  selected  earthquake  list.  Thus  the  data  base  for  this  study 
consists  of  about  100  continental  shallow  local  and  regional 
earthquakes  recorded  by  12  TTVfSSN  stations  in  South  America 
during  the  period  1972*1974,  with  a  range  of  body- wave  magni¬ 
tude  of  4.0  to  6.0. 

The  epicentral  information  was  obtained  from  the  Bulletin  of 
the  International  Seismological  Centre.  In  the  selection  of  the 
data,  only  the  short-period,  vertical-component  microfilm 
records  were  used. 

Figure  5  shows  the  location  of  the  events  used  in  this  study. 
It  also  shows  the  epicenter-station  paths  for  all  the  WWSSN  sta¬ 
tions  and  earthquakes  considered.  As  can  be  seen  in  this  figure, 
most  of  the  paths  are  along  or  across  the  Andes  Mountains  and 
only  a  few  paths  pass  through  the  shield  region.  The  present 
available  data  show  no  usable  earthquakes  in  the  eastern  and 
southern  part  of  the  continent,  even  though  some  parts  of  this 
region  are  covered  by  paths  due  to  earthquakes  which  have 
occurred  in  the  west  and  were  recorded  by  three  stations  (N'AT, 
BDF,  and  LPA)  in  the  east.  Two  big  gaps  (north-northeast  and 
southern  parts  of  continent)  still  exist. 

In  this  study  we  used  first  the  fp  versus  t  method  that  was 
described  earlier  for  different  earthquakes  recorded  at  each  sta¬ 
tion.  In  fact,  to  get  less  scatter  of  the  data  and  a  more  accurate 


value  for  Qt ,  rather  than  finding  an  average  over  a  larger  area, 
the  area  surrounding  each  station  was  divided  into  different 
zones,  and  then  the  predominant  frequency  f?  was  plotted 
against  time  t  after  the  earthquake  origin  for  each  zone.  These 
plots  gave  different  Q,  and  (  values  for  different  directions  about 
the  station.  In  order  to  get  reliable  Qt  values,  it  is  essential  that 
the  observed  predominant  fr  versus  time  t  values  should  cover 
the  sharp  curvature  portion  of  the  master  curves  of  Figure  2> 
Otherwise  there  can  be  large  uncertainties  in  the  Q0  values. 

The  results  of  the  above  process  for  some  stations  are 
presented  and  shown  in  Figures  6  through  8.  In  these  figures, 
circles  denote  the  observed  data  points  and  the  solid  line  is  the 
theoretical  curve  chosen  "by  eye"  to  satisfy  the  data.  The  Q,  and 
f  values  are  also  given  for  each  plot. 

In  each  figure,  two  plots  are  presented  for  each  station, 
which  serve  to  show  the  quality  of  their  fitness  to  the  theoretical 
master  curve.  As  an  example,  in  Figure  3  the  upper  plot  is  for 
the  four  events  north  of  PEL,  whereas  the  lower  plot  is  for  the 
three  events  to  the  south  of  the  station.  Frequency  dependent 
theoretical  curves  for  Q0  *  260  and  <■  =  0.4  (upper  plot)  and  Q0  = 
190,  f  =  0.4  (lower  plot)  are  shown. 

In  Figures  9  through  12  the  events  recorded  by  each  sta.ion 
are  shown  on  the  separate  maps,  along  with  the  Q0  values. 

A  common  feature  of  the  plots  shown  earlier  is  that  they  are 
well  fitted  by  the  theoretical  master  curves,  whereas  for  some 
stations,  like  for  example  TRN  in  Figure  12,  the  scatter  is  cons  id- 
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erable  enough  to  make  it  possible  for  several  curves  with 
different  {  values  to  be  fitted  to  the  data,  resulting  in  different 
and  f  values.  This  is  true  even  if  one  considers  the  shield-type 
paths  only.  This  may  be  due  to  two  possible  factors:  either  the 
areas  designated  as  shields  are  more  heterogeneous  than 
expected,  or  because  the  shield  paths  in  South  America  are  from 
earthquakes  located  in  the  mountains,  the  scattering  in  the  data 
is  caused  by  the  portion  of  the  path  that  crosses  the  mountain 
region. 

To  deal  with  the  problem  of  high  degree  of  scatter  in  the 
data  for  some  stations  like  TRN,  for  which  the  fp  versus  t 
method  fails  to  give  accurate  results,  and  to  check  the  results 
obtained  from  Herrmann's  (1980)  /,  versus  t  method.  Nuttli's  y- 
estimation  method  also  was  used. 

This  method  can  be  simply  applied  to  all  the  earthquakes  by 
taking  one  station  and  one  earthquake  and  measuring  the 
predominant  frequency  and  amplitude  versus  time  after  origin 
time,  starting  from  the  Lg  arrival.  As  explained  earlier,  equa¬ 
tions  (11)  and  (12)  have  been  used  to  obtain  Q0  and  f  for 
different  paths.  For  the  magnitude  correction  factor  (Am*)  in 
equation  (12),  we  used  0. 4-0.6  for  the  western  part.  0  2-0  3  for 
the  northern  part,  and  0.0-0. 1  for  the  eastern  part  of  South 
America.  Figures  13  through  18  present  the  results  obtained  by 
using  Nuttli’s  method. 


Results  of  Q ,  values  obtained  by  Herrmann's  (i960)  and 
Nuttli’s  (1T73)  methods  are  shown  on  Figure  17  and  18.  respec- 


tively.  By  general  analysis  of  these  maps,  one  can  easily  realize 
that  there  are  lateral  variations  in  Qa  across  continental  South 
America.  The  results  obtained  by  application  of  the  two  different 
methods  give  very  much  the  same  Q ,  values  for  the  western  part 
of  the  continent,  whereas  for  the  eastern  part  the  f,  versus  t 
method  gives  higher  values  for  crustal  Qa  than  the  7-estimation 
method.  Figure  19  shows  the  variation  of  the  frequency  depen¬ 
dent  parameter  f  throughout  the  continent,  by  using  coda 
method.  Nuttli(l973)  method  give  higher  f  values  for  eastern 
South  America. 

The  Q4  values  discussed  above  are  obtained  from  long  wave 
paths  (long  travel  times)  which  make  up  the  very  tail  of  the 
seismogram.  To  investigate  whether  these  Q,  values  agree  with 
Q,  values  obtained  from  shorter  wave  paths  (shorter  travel 
times),  a  similar  study  was  done  for  short  wave  paths  with  the 
results  shown  in  Figure  20.  The  important  conclusion  that,  we 
can  get  from  this  comparison  is  that  the  Q,  values  obtained  from 
short  and  long  epicenter  to  station  paths  are  highly  consistent. 
This  map  shows  an  average  Qa  value  of  160  to  220  for  the  western 
part,  a  value  of  550  for  the  southeast  and  an  average  Q0  value  of 
900  for  the  eastern  part.  Figure  21  shows  the  different  zones  of 
Qa  and  frequency  dependent  f  in  South  America. 

Recently  Osagie  and  Mitchell  (1983)  studied  the  attenuation 
of  Qf  in  South  America  at  periods  between  20  and  80s  and  found 
that  there  is  high  attenuation  for  Qfi  in  the  west  of  South  America 
but  lower  in  the  east.  It  is  very  interesting  that,  in  their  study. 


Osagie  and  Mitchell  used  a  simple  crustal  Q  model,  Figure  22,  to 
compute  the  attenuation  and  Q  values  for  1  Hz  Lg  for  South 
America.  Their  numerical  value  is  similar  to  that  of  this  study 
for  the  west  of  South  America.  However,  for  the  east  of  South 
America  they  obtained  a  somewhat  lower  value  for  Qe  than  found 
in  this  study.  The  similarities  of  these  two  studies  are  very 
important,  for  they  provide  more  validation  of  the  results  of  this 
study. 

Finally,  by  using  the  maps  obtained  from  the  fp  versus  t 
method  and  the  7-estimation  methods,  a  contour  map  of  Q„  for 
the  South  American  continental  crust  is  given  in  Figure  23.  This 
figure  represents  one  of  the  important  goals  of  this  study. 

Comparison  of  Qu  Values  with  Tectonics,  Geology  and  Gravity 

The  South  American  continent  can  be  divided  into  three  tec¬ 
tonic  regions,  with  different  origins  and  ages:  the  South  Ameri¬ 
can  Platform,  Patagonian  Platform  and  the  Andean  Cordillera 
Belt.  The  central  and  a  large  part  of  the  eastern  area  of  the  con¬ 
tinent  is  called  the  South  American  Platform,  which  is  the  oldest, 
and  its  basement  is  well  exposed  in  the  three  major  shields 
(Guyana  shield.  Central  Brazil  shield  and  the  Atlantic  shield  in 
the  eastern  part  of  the  continent)  (Figure  24).  Q  values  obtained 
in  this  study  show  that  in  the  Guyana  shield  Q9  varies  from  400  to 
800  and  in  the  Central  Brazilian  shield  from  700  to  1000  (Figures 


17  and  18). 

The  Patagonian  Platform  is  found  entirely  in  Argentina  and 


extends  over  the  broad  continental  margin.  It  is  a  young  plat¬ 
form  (Middle  Paleozoic  onwards).  Q0  values  obtained  for  this 
platform  are  in  the  range  of  420-580  (Figures  17  and  18). 

The  first  two  platforms  are  bounded  by  the  Andean  Cordil¬ 
lera  Belt,  which  is  the  youngest  and  extends  along  the  northern 
boundary  of  the  South  American  Platform.  (Explanatory  note  of 
tectonic  map  of  South  America  by  National  Department  of 
Mineral  Production  of  Brazil  with  the  collaboration  of  UNESCO.) 
Our  Q  map  shows  values  of  115  to  350  for  this  region. 

Correlation  of  Q  Values  -with  Geology  of  the  Continent 

A  geological  map  of  South  America  is  shown  in  Figure  25. 
The  sedimentary  cover  of  the  continent  can  be  compared  with 
the  result  of  our  attenuation  studies  in  South  America.  Gen¬ 
erally,  in  the  active  region  of  the  west  of  South  America  there  is 
no  exact  relation  between  attenuation  and  the  thickness  of  the 
sedimentary  rocks,  whereas  in  the  stable  region  east  of  South 
America,  the  thickness  of  the  sedimentary  rocks  increases  as  Q 
decreases.  This  result  is  similar  to  what  Mitchell  (1983)  observed 
in  the  central  and  eastern  United  States.  Figure  26  shows  the 
variation  of  the  sediment  thicknesses  in  the  eastern  part  of 
South  America. 

Comparison  of  the  Q  Map  with  the  Gravity  Map  of  South  Ame-nca 

The  contoured  Q ,  map  presented  in  this  study  shows  the 
same  pattern  as  the  gravity  map  of  South  America  (Figures  27), 


and  indicates  that  the  low  Qt  values  found  in  the  west  of  South 
America  correlate  with  the  low  gravity  values.  Also,  the  high  Q, 
region  of  the  East  correlates  with  the  high  gravity  anomaly  in  the 
crust.  This  correlation  suggests  that  in  the  region  of  high  tec¬ 
tonic  activity  we  have  high  temperatures  and  high  attenuation, 
and  since  high  temperatures  result  in  rocks  of  low  density,  we 
expect  to  get  relatively  low  Bouguer  gravity  values  in  this 
regions. 

DISCUSSION  AND  CONCLUSIONS 

Previous  Sn  attenuation  studies  in  South  America  indicate 
that,  with  some  exceptions,  there  is  high  attenuation  in  the  crust 
and  upper  mantle  beneath  western  South  America  and  lo’ver 
attenuation  for  the  shield  regions  in  the  central  and  eastern 
parts  of  the  continent  (Molnar  and  Oliver,  1969).  Recent 
attenuation  studies  by  Chinn  at  a L  (1960)  gave  similar  results 
far  Sn  and  Lg  in  South  America.  However,  their  work  >vas 
focused  on  Sn  propagation.  In  this  study,  a  generalized  Qa  map 
for  Lg  waves  of  about  1  Hz  frequency  for  the  South  American 
continent,  as  well  as  contoured  Q ,  maps  of  South  America,  are 
provided  by  using  data  from  12  WWSSN  stations  distributed 
throughout  the  continent,  at  which  over  100  local  and  regional 
earthquakes  with  m*  magnitude  about  4  to  6  were  recorded. 

There  are  two  techniques  which  were  used  for  the  determi¬ 
nation  of  Q:  1)  f9  versus  t  method  proposed  by  Herrmann  (1S80) 
which  is  based  on  a  scattering  model  by  Aki  (1969);  2)  y- 


estimation  method  by  Nuttli  (1973). 

The  Q0  values  for  the  high  frequency  Lg  waves  are  presented 
as  obtained  by  using  both  methods.  The  coda  Qe  values  are  in 
good  agreement  with  the  Qt  values  obtained  from  Lg  amplitudes. 

The  two  different  techniques  should  give  approximately  simi¬ 
lar  results  if  we  have  no  substantial  scatter  in  the  data  and  if 
there  are  sufficient  data.  Although  there  was  some  scattering  in 
the  fp  versus  t  data  for  some  stations,  Q0  values  that  resulted 
from  both  methods  are  highly  consistent  in  the  west  of  South 
America,  a  high  attenuation  region,  but  in  the  east  of  South 
America  the  Q0  values  obtained  by  Herrmann’s  (1980)  method 
are  somewhat  higher  than  those  obtained  using  Nuttli’ s  (1973). 
In  the  western  part  of  the  continent,  except  for  the  region  south 
of  station  PEL,  Q0  data  for  short  epicenter-station  paths  give 
maximum  coverage  of  the  region.  However,  in  eastern  South 
America,  due  to  the  limited  number  of  and  the  epicentral  loca¬ 
tion  of  earthquakes  recorded  by  eastern  stations,  the  Q0  values 
are  averaged  over  larger  areas.  In  addition,  similar  to  the  south¬ 
ern  part  of  the  continent,  the  north-northeastern  section  of  the 
continent  remains  a  gap. 

Data  of  this  study  cover  the  range  of  frequencies  from  0.4  to 
1.4  Hz.  A  frequency  dependence  of  Q  was  investigated  for  this 
range  of  frequencies  by  assuming  Q  =  Q0  (/ /faY- 

Generally,  results  of  this  study  show  there  are  three 
different  zones  of  attenuation  with  different  frequency  depen¬ 
dence  in  South  America,  as  follows: 


1. 


A  high  attenuation  zone  with  large  frequency  dependence,  in 
the  west,  where  Q,  ▼aries  from  130  to  350  and  the 
frequency-dependence  factor  (f)  lies  in  the  range  of  0.4  to 
0.7,  which  carreaponds  to  the  Andean  belt.  This  also  corre¬ 
lates  with  low  gravity  values. 

2.  An  intermediate  attenuation  zone  with  average  frequency 
dependence  in  the  southwest,  central  and  northern  parts  of 
Argentina,  where  Q,  varies  from  420  to  580  and  the  fre¬ 
quency  dependence  factor  (f)  lies  in  the  range  of  0.2  to  0.3, 
which  corresponds  to  the  Patagonian  Platform. 

3.  A  low  attenuation  zone  with  low  frequency  dependence  in  the 
east  of  South  America  with  Qa  of  580  to  980  and  with  a 
frequency-dependence  factor  (0  in  the  range  of  0.0  to  0.2, 
which  corresponds  to  the  South  American  Platform.  This 
also  correlates  with  more  positive  gravity  anomalies. 

Earlier  studies  indicate  that  there  are  higher  for  the  Peru 
region  than  north  of  Chile.  Results  of  this  study  also  agree  with 
higher  Q,  values  for  the  Peru  and  Colombia  regions  than  for  the 
northern  part  of  Chile. 

The  0*  maps  obtained  from  short  and  long  paths  in  this 
study  give  similar  results,  which  is  considered  to  be  proof  that 
the  typical  Lg  coda  is  consistent  with  the  very  late  part  of  the 

coda. 

Recent  results  of  Lg-wave  attenuation  obtained  from  syn¬ 
thetic  seismograms  by  Osagie  and  Mitchell  (1S82)  indicate  simi¬ 
lar  Q,  values  for  the  west  but  lower  Q,  values  for  the  east  than 
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this  study  does. 

The  attenuation  map  for  South  America  is  compared  with 
the  thickness  of  sediment  cover  for  the  continent.  High  attenua¬ 
tion  in  the  west  of  South  America  can  be  correlated  with  6000  to 
15,000  m  of  sediment  thickness.  Generally,  in  the  active  region 
of  the  west  of  South  America  there  is  no  exact  relation  between 
attenuation  and  the  thickness  of  the  sedimentary  rocks,  whereas 
in  the  stable  region  east  of  South  America  Q0  decreases  as  the 
thickness  of  the  sedimentary  rocks  increases.  This  important 
result  is  similar  to  what  Mitchell  (1983)  observed  in  the  central 
and  eastern  United  States. 

As  noted  by  Singh  (1981)  for  North  America,  the  South  Amer¬ 
ican  Q%  values  correlate  with  geological  age,  tectonic  activity  and 
Bouguer  gravity.  In  this  sense  the  contour  map  of  Qa  for  South 
America  (Figure  25)  resembles  the  contoured  Q .  map  of  the 
United  States  obtained  by  Singh  and  Herrmann  (1983). 

The  crustal  Q0  values  given  in  this  study,  as  well  as  the  con¬ 
tour  Qa  map,  might  be  useful  for  studies  of  earthquake  hazard 
and  for  earthquake  engineering  studies  in  South  America. 
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FIGURE  CAPTIONS 

Figure  1.  Cross  section  of  the  Andes  and  the  coastal  region 
(after  James,  1971). 

Figure  2.  Master  curves  for  different  frequency  dependencies  (f) 
for  the  short-period  WWSSN  instruments. 

Figure  3.  /,  vs  t:  Upper  plot  for  four  events  north  of  PEL. 
Lower  plot  for  three  events  south  of  PEL. 

Figure  4.  Time-domain  amplitude  attenuation  of  Lg  waves,  with  7 
as  parameter.  The  units  of  7  are  km~l. 

Figure  5.  Map  showing  the  epicenter-station  paths  for  all  the 
WWSSN  stations  and  earthquakes  considered. 

Figure  8.  /,  vs  t:  Upper  plot  for  four  events  southwest  of  NAT. 
Lower  plot  for  two  events  west  of  NAT. 

Figure  7.  fp  vs  t:  Upper  plot  for  three  events  southeast  of  ANT. 
Lower  plot  for  five  events  southwest  of  ANT. 

Figure  8.  vs  t:  Upper  plot  for  three  events  west  of  TRN 
Lower  plot  for  four  events  southwest  of  TRN 

Figure  9.  Epicenter-station  paths  and  Qa  values  for  earthquakes 
(■)  recorded  at  station  ANT.  Q0  values  are  obtained  by 
using  Herrmann’s  (1980)  method. 

Figure  10.  Epicenter-station  paths  and  Q0  values  for  earthquakes 
(»)  recorded  at  station  PEL.  Q0  values  are  obtained  by 
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using  Herrmann’s  (1980)  method. 

Figure  11.  Epicenter-station  paths  and  G«  values  for  earthquakes 
(•)  recorded  at  station  NAT.  Qe  values  are  obtained  by 
using  Herrmann’s  (1980)  method. 

Figure  12.  Epicenter-station  paths  and  G«  values  for  earthquakes 
(•)  recorded  at  station  TRN.  Q0  values  are  obtained  by 
using  Herrmann’s  (1980)  method. 

Figure  13.  Epicenter-station  paths  and  Qe  values  for  earthquakes 
(•)  recorded  at  station  ANT.  Qq  values  are  obtained  by 
using  Nuttli’s  (1973)  method.  (Refer  to  Figure  9  for  the 
coda  Q  estimates.) 

Figure  14.  Epicenter-station  paths  and  Q0  values  for  earthquakes 
(•)  recorded  at  station  NAT.  Q0  values  are  obtained  by 
using  Nuttli’s  (1973)  method.  (Refer  to  Figure  11  for 
the  coda  Q  estimates.) 

Figure  15.  Epicenter-station  paths  and  Go  values  for  earthquakes 
(•)  recorded  at  station  PEL.  Go  values  are  obtained  by- 
using  Nuttli’s  (1973)  method.  (Refer  to  Figure  10  for 
the  coda  Q  estimates.) 

Figure  16.  Epicenter-station  paths  and  Go  values  for  earthquakes 
(«)  recorded  at  station  TRN.  Go  values  are  obtained  by 
using  Nuttli's  (1973)  methoci.  (Refer  to  Figure  12  for 
the  coda  Q  estimates.) 


Figure  17.  Average  crustal  Q ,  values  for  the  7TWSSN  stations  by 
using  Herrmann's  (1980)  method. 

Figure  18.  Average  crustal  Q ,  values  from  the  WWSSN  stations  by 
using  Nuttli’s  (1973)  method. 

Figure  19.  Distribution  of  {  values  throughout  the  South  Ameri¬ 
can  continent,  by  using  coda  method. 

Figure  20.  Average  crustal  Q,  values  from  the  WWSSN  stations 
obtained  from  shorter  wave  paths  using  Herrmann’s 
(1980)  method. 

Figure  21.  Map  showing  different  zones  of  Q  and  frequency 
dependence  of  f  in  South  America. 

Figure  22.  Models  of  anelasticity  for  eastern  and  western  South 
America  (from  Osagie  and  Mitchell,  1983). 

Figure  23.  Contour  map  of  crustal  Q,  for  the  entire  South  Ameri¬ 
can  continent. 

Figure  24.  Tectonic  map  of  South  America  compiled  by  Kent  C 
Condie  (1982). 

Figure  25.  Geologic  map  of  South  America  compiled  by  Gus  H 
Goudarzi  (1977). 

Figure  28.  Map  showing  the  variations  of  the  sediment 
thicknesses  in  eastern  South  America. 


Figure  27.  Bouguer  gravity  map  of  South  America  showing  the 
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Effeot  of  Crustal  Velooity  Structure  and 
Crustal  Q  Values  on  the  Amplitudes  and  Wave 
Forms  of  Lg 
by 

B.J.  Mitchell  and  H.J.  Hwang 

Introduction 

The  amplitudes  and  wave  forms  of  the  phase  Lg  are  influenced  by 
several  factors,  including  velocity  and  Q  structure  of  both  crustal  and 
sedimentary  rock.  Variations  in  these  factors  with  depth,  as  well  as 
lateral  complexities  along  the  path  of  travel  are  expected  to  play  a 
role  in  the  wave  farms  recorded  at  seismic  stations.  Our  work  over  the 
period  of  this  contract  has  been  concerned  with  studying  those  factors 
which  can  be  addressed  using  plane- layered  models  of  the  crust.  Thus 
our  results  either  ignore  the  effect  of  lateral  complexities  or  consider 
them  to  be  approximately  modelled  by  layered  models  over  their  paths  of 
travel. 

In  order  to  minimize  the  effect  of  lateral  complexities,  we  have 
used  two  approaches.  First,  we  have,  as  far  as  possible,  used  rela¬ 
tively  short  paths  which  avoid  major  lateral  boundaries,  such  as  con¬ 
tinental  margins  or  major  boundaries  between  tectonic  provinces. 
Second,  we  have  used  the  early  portions  of  observed  Lg  wave  forms  in  the 
data  analysis  portion  of  our  study.  Thus  we  have  used  waves  which  are 
more  coherent,  and  presumably  less  scattered,  than  the  later  portions  of 
the  Lg  coda. 

In  the  eastern  United  States,  velocity  and  Q  models  have  been  found 


which  adequately  predict  observed  surface  wave  phase  velocities,  group 
velocities  and  attenuation  at  both  intermediate  and  short  periods. 
Models  have  also  been  found  which  satisfy  similar  data  for  other 
regions,  suoh  as  the  Basin- and-Range  province,  although  in  the  case  of 
Q,  the  models  are  less  reliably  determined.  We  have  studied  various 
factors  which  affect  the  attenuation  of  higher  mode  Rayleigh  waves  (or 
Lg)  in  various  period  ranges.  We  will  first  discuss  the  effect  of  the 
velocity  model  on  amplitudes  and  attenuation  in  both  the  eastern  United 
States  and  Basin- and-Range.  Second,  we  will  attempt  to  see  whether  or 
not  the  observed  regional  variation  of  1-Hz  Lg  in  both  eastern  and 
we:  tern  United  States  can  be  explained  simply  by  the  variations  in 
thicknesses  of  sediments  there.  Third,  we  will  test  two  Q  models  of  the 
Basln-and-Range  province  to  see  if  the  rapid  attenuation  known  to  occur 
there  can  be  explained  by  low-Q  sediments  over  a  Q-model  similar  to  that 
for  the  stable  regions  in  the  eastern  United  States,  or  if  it  requires 
low  Q  values  in  the  deeper  crust. 

Effect  of  Veloaity  Structure  on  Attenuation  and  Yield 

Four  velocity  models  derived  from  data  in  the  central  and  eastern 
United  States  were  used  to  compute  short-period  Lg  synthetics.  These 
are  a  simplified  version  of  a  model  by  Mitchell  and  Herrmann  (1979),  as 
well  as  models  developed  by  McEvllly  (1964),  Herrmann  (1979),  and  Taylor 
(1980).  These  are  shown  in  Table  1.  The  simplified  eastern  United 
States  model  in  that  table  (Mitchell,  1982)  is  composed  of  six  layers, 
and  was  used,  rather  than  the  multi-layer  model  of  Mitchell  and  Herrmann 
(1979),  in  order  to  save  computation  time. 

All  of  the  velocity  models  were  combined  with  a  Q  model  in  which 
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coopressional  wave  Q  ( Q , )  la  2000  and  shear  wave  Q  (Q„)  is  1000  at  all 

d  P 

depths.  Short-period  synthetics  for  these  models  were  calculated  using 
programs  of  Wang  (1980).  Our  results  indicated  that  only  those  models 
which  include  a  surface-layer  of  low-velocity  sediments  produce  realis¬ 
tic  looking  seismograms.  Those  models  without  sediments,  and  comprised 
of  thick  crustal  layers,  contain  spike-like  arrivals  which  are  not 
observed  on  real  seismograms.  Examples  for  two  models  appear  in  Figures 
la  and  1b.  The  variations  in  velocity  structure  can  produce  somewhat 
different  amplitudes  at  a  given  distanoe  from  the  source  and  can  lead  to 
small  differences  in  magnitude,  and  consequently  in  yield,  for  this  case 
in  which  the  source  is  located  at  a  depth  of  5  km. 


Synthetics  were  next  computed  for  a  velocity  model  of  the  Basln- 

and-Kange  province  (Baohe  et  ql.,  1978)  and  these  were  compared  with 

those  computed  for  the  simplified  eastern  United  States  model.  These 

models  appear  in  Figure  2.  Identical  Q  models  were  used  for  both 

regions,  Q_  being  set  at  300  through  the  upper  18  km  of  the  crust  and 

P 

1000  at  greater  depths.  Q.  was  assumed  to  be  twice  as  large  as  Q_  at 

«  p 

all  depths. 


Table  2  gives  magnitude  differences  produced  by  differences  in 
velocity  structure  between  the  simplified  EUS  and  Basin-and-Range  model 
for  a  source  at  a  depth  of  0.5  km.  The  source  corresponds  to  an  explo¬ 
sion  in  the  low- velocity  sediment  layer  of  each  model.  Although  the 
magnitude  differences  for  1-Hz  Lg  are  as  large  as  0.3  at  a  distance  of 
750  km.  they  become  less  consequential  at  greater  distances.  Table  2 
also  lists  differences  in  m.  derived  from  spectral  amplitudes  of  larger 

D 

period  values.  At  periods  of  8  and  10  s,  these  differences  are  very 


small,  although  they  are  larger  at  short  periods. 

Effect  of  Sediments  on  Q.  in  the  Eastern  United  States 

Lg 

Q  values  obtained  from  Lg  amplitudes  (Espinosa,  1981)  and  from  the 
dispersion  of  the  Lg  coda  (Singh  and  Herrmann,  1 983)  exhibit  a  similar 
regional  variation  across  the  United  States.  The  map  of  Singh  and 
Herrmann  shows  the  highest  Q  values  in  the  north  and  central  United 
States,  with  gradual  decreases  extending  both  eastward  and  westward  from 
there  into  the  mid-Atlantic  states  and  Great  Plains.  More  rapid 
decreases  oocur  along  the  Atlantic  and  Gulf  coasts. 

The  distribution  of  sediments  everywhere  east  of  the  Rocky  Moun¬ 
tains  (e.g.  Sloss  fit  al. ,  I960),  especially  that  of  young  sandstones  and 
shales,  bears  a  resemblance  to  the  areal  distribution  of  the  G  values 
from  Lg  coda  waves. 


In  an  attempt  to  see  if  the  regional  variation  of  Q  for  Lg  waves 
could  be  explained  simply  by  the  effect  of  sediment  distribution,  an 
attempt  was  made  to  model  those  regions,  using  the  crustal  Q  model  of 
Mitchell  (1980)  overlain  by  various  thicknesses  of  sediments,  depending 
upon  geographic  region.  It  was  assumed  that  crustal  Q  was  frequency- 
dependent  whereas  the  Q  values  of  the  sediments  were  not.  It  was  also 
assumed  that  Q  values  of  limestones  and  dolomites  were  similar  to  those 
of  crustal  rock  and  that  in  the  low-Q  sediments,  sandstone  and  shales. 


Figure  3  (from  Sloss  fit  al. ,  I960)  illustrates  the  distribution  of 
sediments  in  the  United  States  for  the  Cretaceous  period.  Shales  and 
sandstones  for  that  period  thicken  to  the  west  and  along  the  east  and 
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Gulf  coasts  as  compared  with  the  north-central  United  States.  Sediments 
from  other  periods  exhibit  other  patterns.  Combining  these  thicknesses 
with  sediments  of  various  other  periods,  we  constructed  models 
corresponding  to  regions  in  western  Pennsylvania,  Nebraska,  and  southern 
Louisiana.  Values  of  Q  were  assigned  to  those  sediments  using  the  Q 
versus  pressure  curve  of  Winkler  and  Nur  (1979)  for  partially  saturated 
sandstone.  Table  3  presents  simplified  Q  models  of  sediments  for  the 
various  regions.  These  were  combined  with  the  frequency-dependent  Q 
models  of  Mitchell  (1980)  for  the  eastern  United  States.  Table  4  com¬ 
pares  Q  values  for  Lg  computed  for  those  models,  as  well  as  for  a  model 
of  the  high-Q  region  in  the  central  United  States,  and  compares  them 
with  values  observed  by  Singh  and  Herrmann  (1983). 

Our  results  Indicate  that  it  is  very  easy  to  obtain  Q  values  for  Lg 
waves  which  are  as  low  or  lower  than  those  which  have  been  observed, 
simply  by  using  realistic  Q  values  for  sediments.  In  fact,  our  first 
results  yielded  Q  values  for  Lg  which  were  too  low  in  all  cases.  These 
Q  differences  produced  by  regional  variations  of  low-Q  sediments  can 
produce  very  large  differences  in  m^.  Table  5  shows  that  these  can  be 
as  large  as  0.7  in  severe  cases. 

These  results  suggest  that  Lg  Q  in  stable  regions  can  be  estimated 
if  we  have  a  knowledge  of  near-surface  geology. 

Effect  of  Sediments  on  Q.  in  the 

Lg 

Western  United  States 

Results  of  computations  in  the  previous  section  indicate  that  the 
rather  large  regional  variation  of  Q  observed  in  the  eastern  and  central 
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United  States  can  be  explained  as  being  due  In  large  part  to  variations 
In  the  thickness  of  surflcial  sediments.  Variations  of  Lg  Q  values  also 
occur  in  the  western  Onlted  States,  although  to  a  smaller  extent.  It 
will  be  important  to  know  whether  or  not  those  variations  are  also  pro¬ 
duced  by  varying  accumulations  of  sediment. 

He  first  found  a  crustal  model  which  produced  Lg  wave  forms  which 
are  similar  to  those  observed  for  paths  across  the  Basin-and-Range  pro¬ 
vince.  The  model  found  to  be  most  satisfactory  lnoluded  a  low-Q  upper 
crust  (Q^  a  Q  a  150)  and  a  hlgher-Q  lower  crust  (Q^  =  a  500).  This 
Q  model  was  combined  with  the  velocity  model  of  Chen  and  Mitchell  ( 1984) 
for  that  region.  The  Lg  synthetlo  seismograms  produced  by  that  model 
are  shown  in  Figure  4  for  4  distances. 

Varying  thicknesses  of  sediments  (1,  2,  and  4.  km)  were  placed  above 
that  model  as  shown  in  Table  6  to  see  how  that  variation  might  affect  Lg 
Q  values.  As  in  the  previous  section.  Lg  Q  values  were  taken  to 
increase  with  pressure  in  accordance  with  laboratory  measurements. 

Theoretically  expected  Q  values  for  Lg  are  shown  beneath  each 
model.  The  maximum  range  is  between  190  and  230.  These  results  indi¬ 
cate  that,  in  contrast  to  stable  hlgh-Q  regions,  such  as  the  eastern 
United  States,  Lg  Q  values  for  regions  with  low-crustal  Q  values  are  not 
greatly  affected  by  varying  thickness  of  sediment,  if  losses  due  to 
scattering  are  ignored. 

Q  Structure  in  the  Basin-and-Range  Province 

Mitohell  (1975),  in  inverting  intermediate-period  attenuation  data 
for  the  western  United  States,  found  that  Q_  values  in  the  upper  crust 

P 


in  the  vestern  United  States  were  much  lower  than  those  in  the  eastern 


United  States.  Cheng  and  Mitchell  (1981),  using  a  new  multi-mode  method 
and  paths  nearly  totally  confined  to  the  Basin-and-Range  province  found 


very  low  Q_  values  (approximately  85)  in  the  upper  crust.  Later  work 
P 


has,  however,  suggested  that  it  may  be  possible  to  explain  the  fall-off 
of  amplitudes  at  intermediate  periods  by  thick  layers  of  low-Q  sediment 
overlying  hlgh-Q  crust. 


Two  types  of  models  were  tested  using  both  long-period  and  short- 
period  seismograms.  Long-period  seismograms  computed  for  the  Bache  fit 
al.  ( 1978)  velocity  model  and  the  two  Q  models  are  nearly  identical  in 
form,  so  single  station  data  at  long  periods  cannot  be  used  to  distin¬ 
guish  between  the  two  models. 


The  fall-off  of  1-Hz  Lg  amplitudes  was  then  determined  from  short- 
period  synthetic  seismograms  predicted  by  the  two  models.  The  Q  values 
for  Lg  in  these  cases  are  550  for  the  model  with  the  high-Q  crust  and 
low-Q  sediments  and  290  for  the  model  with  the  low-Q  crust.  Thus  it 
appears  that  the  low  Q  values  of  Lg  in  the  Basin-and-Range  province  can¬ 
not  be  due  to  low-Q  sediments,  but  reflect  low  Q  values  in  the  upper 

crust.  Moreover,  a  model  in  which  Q_  varies  with  frequency  is  not 

P 

required  to  explain  both  Intermediate-  and  short-period  data. 
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Figure  la.  Vertical  component  synthetic  seismograms  at  4  distances 
for  the  Simplified  EUS  model.  Qa  -  2000  and  =  1000  at  all 
depths. 


115 


JO 


oO 


90 


ISO 


1  50  1  SO 

T  *  300.00 


2i  0 


240 


270 


JCO 


3J0 


•  40  150  180  ~2\0  240  sTo  300 

T-4C0.00 


3  JO 


Figure  lb.  Vertical  component  synthetic  seismograms  at  4  distances 
for  the  velocity  model  of  McEvilly  (1964).  Q  =  2000  and 
Q$  ■  1000  at  all  depths.  a 
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One  of  the  first  observations  to  be  made  when  regional  studies  of 
attenuation  began  in  a  systematic  way  in  the  mid-70's  was  that  shear 
wave  Q  values  (Q_)  in  the  crust  obtained  from  intermediate-period  sur- 

P 

face  waves  can  vary  by  large  amounts  from  region  to  region  (Mitchell, 
1975;  Cheng  and  Mitchell,  1981).  Q  values  for  1-Hz  Lg  waves  were  also 
found  to  vary  regionally  in  much  the  same  way  that  Q_  from  surface  waves 

P 

were  found  to  vary.  Studies  such  as  that  of  Singh  and  Herrmann  (1983) 
presented  maps  of  1-Hz  Q  values  which  showed  low  values  of  Q  for  Lg 


waves  in  the  same  region  where  low  values  of  Q_  in  the  upper  crust  had 

P 


been  found. 


Q,  has  also  been  found  to  be  frequency-dependent  in  the  upper  crust 

P 

(Mitchell.  1980).  The  degree  of  that  frequency  dependence  varies 


regionally,  being  higher  in  regions  of  high  Q_  and  lower  in  regions  of 

P 

low  Q_  in  the  United  States  (Mitchell,  1981). 

P 


The  purpose  of  the  present  report  is  to  discuss  preliminary 
attenuation  data  obtained  at  intermediate  periods  from  South  America  and 
from  India.  Both  regions  have  large  areas  which  are  tectonically  stable 
as  well  as  large  areas  which  have  undergone  extensive  tectonic  activity. 
It  should  therefore  be  possible  to  verify  whether  our  results  for  stable 
and  aotive  regions  of  North  Amerioa  are  also  applicable  to  other  con¬ 


tinents 
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Data  Acquisition  and  Processing 

Seismograms  which  included  well-developed  fundamental-mode  Rayleigh 
waves  were  obtained  for  several  earthquakes  in  both  South  America  and 
India.  The  events  were  selected  such  that  they  were  situated  on  a  com¬ 
mon  great  circle  path  between  two  seismograph  stations. 

For  the  selection  of  events,  a  computer  program  was  used  to  search 
the  tape  of  OSGS  PDE  epicenters  for  events  located  within  4  degrees  of 
arc  of  the  azimuth  of  a  line  Joining  two  seismograph  stations.  The  pro¬ 
gram  obtains  a  set  of  earthquakes  for  each  of  the  station  pairs  used.  A 
manual  selection  was  then  made  from  this  set  of  events  on  the  basis  of 
magnitude,  distance  to  the  station,  and  depth.  The  data  used  were  vert¬ 
ical  component  seismograms  recorded  by  instruments  of  the  World-Wide 
Standard  Seismograph  Network  (WWSSN). 

Spectral  amplitudes  and  group  velocities  of  the  fundamental-mode 
Rayleigh  wave  were  obtained  by  the  multiple-filter  technique  (Dziewonaki 
Alt.  gJL. ,  1969).  The  recorded  seismograms  were  digitized  at  irregular 
intervals  and  a  sampling  interval  of  1.0  sec  was  then  obtained  by  linear 
interpolation  of  the  digitized  points.  Digitization  started  1  or  2 
minutes  before  the  arrival  of  the  surface  waves  and  extended  to  times 
corresponding  to  group  velocities  of  about  2.2  km/ sec.  The  length  of 
the  digitized  record  depends  on  the  length  and  character  of  the 
dispersed  wave,  but  most  often  it  is  8-10  minutes  long. 

Attenuation  at  Intermediate  Periods 

Attenuation  coefficients  were  computed  from  amplitude  spectra  using 
the  two-station  method  (Ben-Menahem,  1965).  An  advantage  of  the  two- 


station  method  is  that  it  does  not  require  a  knowledge  of  the  source  or 
the  velocity  structure  along  the  path  of  propagation.  Several  indivi¬ 
dual  determinations  along  the  same  path,  however,  can  show  large  fluc¬ 
tuations  due  to  multlpathing.  The  attenuation  coefficients  are  obtained 
using  the  equation 

(1) 

where  A  refers  to  speotral  amplitude.  Z  to  instrument  response,  r  to 
distance  in  km,  and  £  to  distance  in  degrees.  The  subscripts  1  and  2 
refer  to  values  at  the  near  and  more  distance  stations,  respectively. 

Preliminary  Rayleigh  wave  attenuation  coefficients  determined  for 
the  eastern  and  western  portions  of  South  America  appear  in  Figure  1. 
The  scatter  in  the  data  is  reflected  in  standard  deviation  bars  which 
accompany  the  mean  values  at  each  period.  It  is  quite  apparent  that  the 
attenuation  coefficients  for  western  South  America  are  substantially 
higher  than  those  for  eastern  South  America,  especially  at  shorter 
periods  between  20  and  40  seconds.  The  values  for  eastern  South  America 
are  similar  to.  and  perhaps  slightly  lower  than,  those  previously  deter¬ 
mined  for  eastern  North  America.  The  values  for  western  South  America 
are  similar  to  those  obtained  in  tectonically  active  regions  of  North 
America  (Mitchell  atgl.,  1977)  or  Eurasia  (Yacoub  and  Mitchell.  1977). 
Similar  attenuation  coefficient  values  have  been  obtained  for  the  stable 
and  tectonically  active  regions  of  India  from  the  portion  of  the  avail¬ 
able  data  from  that  region  which  has  been  analyzed  up  to  now. 
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Inversion  of  Q; 


P 


The  equations  of  Anderson  et  gl*  (1965)  were  used  to  invert  the 
attenuation  ooeffioient  data,  assuming  that  Q_  is  independent  of  fre- 

P 

quency.  In  order  to  utilize  those  equations,  a  velocity  aodel  is 
required  for  the  region  of  study.  Phase  velocities  were  obtained  using 
the  two  station  method.  Across  the  eastern  region  of  South  America, 
they  are  very  close  to  the  phase  velocities  of  the  Canadian  shield 
(Brune  and  Dorman.  1963)*  The  phase  velocities  determined  for  the 
western  portion  of  South  America  agree  closely  with  values  determined 
for  that  region  by  James  (1971).  Because  of  the  agreement  of  phase 
velocities  observed  in  the  present  research  with  those  obtained  from 
earlier  models,  we  use  those  veloolty  models  to  compute  partial  deriva¬ 
tives  to  be  used  in  the  inversion  process  for  Q_. 

P 

Models  of  Q_  were  obtained  for  both  eastern  and  western  South  Amer- 

P 

lea  (Figure  2).  The  Q_  values  for  western  South  America  appear  to  be 

P 

lower  at  all  depths  than  those  for  eastern  South  America,  although  at 
greater  depths  in  the  crust  it  is  not  possible  to  resolve  the  differ¬ 
ences  between  the  two  regions.  The  data  for  western  South  America 
appear  to  require  a  low  Q  zone  in  the  upper  mantle  whereas  the  data  for 
eastern  South  America  do  not.  Because  of  the  similarity  of  the  attenua¬ 
tion  data  for  India,  it  is  also  expected  that  crustal  Qn  values  will  be 

P 

lower  in  tectonloally  active  regions  than  in  stable  regions  of  that 
sub-continent. 


Attenuation  of  1-Hz  Lg 


The  models  shown  in  Figure  2  were  used  to  compute  the  attenuation 
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and  Q  values  for  1-Hz  Lg  and  to  oompare  them  with  previously  determined 
observed  values  (Raoof  and  Nuttli.  1983) »  Programs  written  by  C.Y. 
Vang  and  R.B.  Herrmann  (Vang,  1981)  were  used  to  compute  synthetic 
seismograms,  corresponding  to  short-period  Instruments  of  the  Vorld-Vide 
Standard  Seismograph  Network  (WSSN),  and  to  measure  the  amplitude 
fall-off  with  diatanoe  to  determine  attenuation  coefficients  and  Q 
values.  This  method  allows  us  to  determine  Q  values  for  Lg  predicted  by 
various  models  in  muoh  the  same  way  that  Q  is  determined  from  real 
seismograms. 


Synthetic  seismograms  of  1-Hz  Lg  waves  were  computed  for  the  models 
of  both  eastern  and  western  South  America  (Figures  3  and  4).  The 
seismograms  for  western  South  America  not  only  fall  off  more  rapidly 
with  distance  than  those  for  eastern  South  America,  but  also  have  lower 
predominant  frequencies,  especially  at  greater  distances.  The  models 
predict  values  of  Q  for  Lg  recorded  on  short-period  WSSN  instruments  to 
be  480  and  290.  respectively,  for  eastern  and  western  South  America. 


The  Q  value  of  Lg  predicted  by  the  model  of  eastern  South  America 
(480)  is  substantially  lower  than  those  observed  by  Raoof  and  Nuttli 
(1983).  This  result  implies  that  Q_  in  the  crust  beneath  eastern  South 

P 

America  is  substantially  higher  at  high  frequencies  (near  1-Hz)  than  at 
lower  frequencies.  This  result  is  similar  to  that  of  Mitchell  (1981) 
for  eastern  North  America.  By  contrast,  the  Q  value  of  Lg  predicted  by 
the  model  of  western  South  America  is  similar  to  that  observed  by  Raoof 
and  Nuttli  (1983).  Qa  in  the  crust  beneath  western  South  America. 

P 

therefore,  does  not  vary  significantly  with  frequency  over  the  frequency 


range  considered  in  this  report.  These  results  are  consistent  with  ear- 


2 


Her  results  In  North  America  (Mitchell.  1981),  where  values  in  the 
crust  of  the  tectonically  aotive  region  of  western  North  America  do  not 
vary  greatly  with  frequency  whereas  they  do  in  the  stable  craton  of 
eastern  North  America.  Preliminary  results  for  stable  and  tectonically 
aotive  regions  of  India  also  seem  to  be  consistent  with  that  pattern. 
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Figure  3.  Synthetic  seismograms  for  paths  across  eastern  South  America  for 
the  phase  Ig  as  recorded  on  short-period  WWSSN  seismographs. 
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Figure  4.  Synthetic  seismograms  for  paths  across  western  South  America  for 
the  phase  Lg  as  recorded  on  short-period  WWSSN  seismographs. 


